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TECHNICAL MEMORANDUM 


THEORY OF OSTWALD RIPENING IN A TWO-COMPONENT SYSTEM - 
CENTER DIRECTOR’S DISCRETIONARY FUND FINAL REPORT 


I. INTRODUCTION 

1.1 Physical Picture of Ostwald Ripening 

L* 

* Imagine a solution containing two chemical components. The solution is 

supersaturated with solute, and nuclei of a second phase have already started 
to grow. The supersaturated solution shall be called the parent phase, and 
the nuclei shall be termed the second phase. The nuclei shall be assumed 
f large enough to have well-defined surfaces so that macroscopic theories such 

as thermodynamics can be applied. Such nuclei shall be called grains. Of 
interest is the grain size distribution function and its development over 
time. The grain size distribution evolves so as to reduce the surface area of 
contact between the two phases, which also serves to reduce the interfacial 
free energy. The dynamics of this process is termed Ostwald ripening [1]. 
Here, it shall be assumed that the transfer of material between the two phases 
is governed by Fick's laws of diffusion. The two phases may be any pairing of 
solid, liquid, or gas, except both cannot be gases, since gases are miscible 
in all proportions. 

During Ostwald ripening, the solute super saturation in the parent phase 
diminishes with time tending asymptotically to zero. At any given time, the 
value of the super saturation determines a critical grain size. Grains larger 
than this tend to grow, while grains smaller tend to dissolve. The parent 
phase is the medium through which the grains communicate by way of mass trans- 
fer. Associated with each grain are local concentration gradients of solute 
and solvent. For example, a grain which is growing by concentrating solute 
and rejecting solvent has about it a positive solute gradient and a negative 
solvent gradient. 

It is a main objective of this report to treat quantitatively the 
situation where diffusion of both solvent and solute occur in the parent 
phase. Associated with these diffusion fluxes are bulk transfers of material 
termed advections. These shall be treated by referring the diffusion to a 
coordinate frame in the parent phase which moves with respect to the labora- 
tory [2] . Our result will permit the extension of the current theory of 
Ostwald ripening [3-5] to the case where both chemical components occur in 
** both phases. Phase separating systems possessing a miscibility gap are an 

important example. 

As the total volume occupied by the grains increases, the grains come 
closer together, and the associated concentration gradients begin to overlap. 
Overlapping concentration gradients mean interacting local rates of growth and 
dissolution. The rate of growth of an individual grain is no longer the same 
as it would be if it were in isolation; instead, many-body effects come into 
play. 



It is a second main objective of this report to incorporate these many- 
body effects into the theory of Ostwald ripening. The work of Marqusee and 
Ross [6] will be followed to accomplish this. An extended, pedagogical de- 
velopment of their original ideas shall be given. 

Our program consists essentially of four tasks: 

( 1 ) Derive the Gibbs-Kelvin equation relating the solubility of a grain 
to its size (see Section 1.2). 

(2) Taking into account growth and dissolution, develop a continuity 
equation for the time evolution of the grain size distribution function (see 
Section 1.3). 

(3) At any given time the total masses of solute and solvent are di- 
vided between the parent phase and the second phase. This division must 
satisfy a conservation equation, derived in Section 1.4. 

(4) In the continuity equation, there occurs a function depending upon 
the grain size and the time, which specifies the rate of growth of an indi- 
vidual grain. This growth velocity function, which is derived in Section 1.5, 
depends on a carefully defined diffusion coefficient. All of Chapter II is 
occupied with the problem of generalizing this function to take into account 
the effect of overlapping concentration gradients. This technically involved 
subject will be treated within the framework of the classical many-body theory 
of statistical mechanics. 

In summary. Chapters I and II are devoted to the program above. The 
work of combining the results of this program into a single differential 
equation is presented in Sections 3. 1-3.4 of Chapter III. In Section 3.5, 
that equation is integrated, and the functional form of the grain size 
distribution function is found. Parameters associated with that distribution 
are determined as a function of the volume fraction of the system in Sections 
3. 6-3. 7. In Section 3.8, results are presented graphically and by equations 
in a form ready for comparison with experiment. The system 

succinonitrile/water is analyzed as an example. At the end of the report are 
five appendices explaining the mathematical techniques used in Chapters I-III 
and deriving certain technical details associated with the solution of the 
continuity equation for the grain size distribution function. 

1 • 2 Solubility and the Gibbs-Kelvin Equation 

A molecule in the interior of a precipitate grain is bound by forces 
which are characteristic of the bulk precipitate phase. A molecule in the 
surface of the grain, however, has some of its atomic bonding contacts satis- 
fied by molecules in the parent phase which are adjacent to it. These con- 
tacts provide less stabilizing (negative) energy for the surface molecule than 
the contacts it would enjoy if it were in the bulk of the grain. Hence, in an 
algebraic sense the chemical potential of a moleucle in the interior of the 
grain is less than a molecule in the surface. 
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A 




Suppose now that the parent phase is an ideal solution composed of 
components A and B with concentrations ( molecule s/cm^ ) C A and Cg, respective- 
ly 

are given by 


The corresponding chemical potentials are y and y , respectively, which 

A B 


"A - "A + Mln C A 


( 1 . 2 . 1 ) 


»B * "B * kT1 " °B 


( 1 . 2 . 2 ) 


where y and y are the chemical potentials defined in the standard states C A 
B B A 

= 1 and Cg = 1 . Here k is Boltzmanns constant and T is the absolute tempera- 
ture. 


Consider a grain which is part of the second phase that has precipitated 

V I 

fran the parent phase. Let there be n A and n B molecules of A and B in this 

grain and let its surface area be S’ and its volume be V 1 . The composition of 

the second phase is, of course, different than that of the parent phase. The 

chemical potentials of A and B molecules in the bulk of the grain are y ° 

« A 

1 O 

and y B , respectively. The grain forms a surface with the parent phase with 
surface energy given by y. The Gibbs free energy of component A in the grain 
is 


G A = VV° + ys’ (1.2.3) 

and of component B 

°B = Vb° + yS ' • (1.2.4) 

To obtain chemical potentials for the grain as a whole, we must dif- 

« V 

ferentiate Eqs. (1.2.3) and (1.2.4) with respect to n^ and n^, 
respectively. In so doing, we should hold in mind that S* depends upon the 

i t 

number of molecules, n^ and n R , in the grain. In other words, the larger the 
grain, the larger its surface area. The chemical potential of A* in the grain 
is 


3 



3G 

1 

1 

A _ 

' o , dS 

3V 

i 

U A + Y — T 

# i 

8n A 

dV 

8 "a 


(1.2.5) 


I « 

where to compute 3S /9n we have used the chain rule with V* as the inter- 

A 

mediate variable. 

The derivative dS'/dV depends on the geometric shape of the grain. At 

2 

this point we specialize our analysis to spheres for which S 1 = 4 ttR and V' = 
3 

(4tt/ 3)R where R is the radius. The derivative is 


i 


i 

dS _ 87rRdR _ 2 

1 2 

dV 4 ttR dR R 


* 

( 1 . 2 . 6 ) 


The partial molecular volume of component A in the grain is 


« 




(1.2.7) 


Substituting Eqs. (1.2.6) and (1.2.7) into Eq. (1.2.5), we obtain 


I 



At equilibrium = y^. This determines the dependence on C A on R. 
Equating Eqs. (1.2.1) and (1.2.8), we obtain 


kTln C a (R) = -(y° 




2 J\ 

R 


(1.2.9) 


V 


✓ 


In the limit R-h», C (R) + C (») , and according to Eq. (1.2.9) 

A A 

kTln C A (~) = -(u A - y A °) . (1.2.10) 

Equation (1.2.10) gives the solubility, C A (°») , of a flat surface (infinite 
radius of curvature) and is the usual equilibrium expression for the concen- 


4 


tration of a component distributed between two phases. We substitute Eq 
(1.2.10) into Eq. (1.2.9) to obtain 


* 


* 


w 


C (R) = C (<*>) exp( 2yv' /kTR) . 

A A A 


( 1 . 2 . 11 ) 


Equation (1.2.11) is called the Gibbs-Kelvin equation and shows how the equi- 
librium solubility of a component outside a curved surface depends upon 
radius . 

Starting from Eqs. (1.2.2) and (1.2.4), we may derive by analogy 

C (R) = C (co) exp(2^v'/kTR) . (1.2.12) 

5 B B 

1.3 Continuity Equation 

Consider a system at time t consisting of a large number of grains. Let 
n(R,t)dR be the number of grains per unit volume with radii lying between R 
and R + dR. As the grains grow in size, their number must decrease because 
the total mass in the system must be conserved. Suppose that the time rate of 
change of the radius of any grain is known and is expressed by the equation 


dR(t) 

dt 


V(R(t) ,t) 


(1.3.1) 


where the function V(R(t),t) contains all of the information specifying the 
growth dynamics. Suppose further that at a fixed time t, we are able to 
determine the grain densities n(R(t),t) and n(R(t) + dR,t) at radii R(t) and 
R(t) + dR, respectively. The flux of grains leaving R(t) is n(R(t),t) 
V(R(t),t) while the flux leaving R(t) + dR is n(R(t) + dR,t)V(R(t) + dR,t) • 
The net rate of buildup of grains between R(t) and R(t) + dR is (3n(R(t),t)/ 
3t)dR. To conserve numbers, we must have 


9n(R(t),t) 

3t 


dR = n(R(t) , t) V(R( t) ,t) 


- n(R(t) + dR,t)V(R(t) + dR, t) 


(1.3.2) 


We expand the last term on the right-hand side of Eq. (1.3.2) in a Taylor's 
series about R(t), 

5 


n(R( t) + dR,t)V(R(t) + dR, t) = n(R(t) ,t)V(R(t) ,t) 

. 9 ( n(R( t) ,t)V(R(t) ,t)) 

+ 8R dR * (1 * 

Substitution of Eq. (1.3.3) into Eq. (1.3.2) and cancellation of dR yields 


3n(R(t),t) . 9( n(R( t) , t)V(R( t) , t) ) „ .. . .. 

3t + 9R ° * d.3.4) 

Equation (1.3.4) is called the continuity equation for the grain size density 
n(R(t) ,t) . 

A model for the "velocity,” V(R(t) ,t) , will be developed in Section 1.5. 
1.4 Conservation of Mass 

At any time t, the total mass of any component (say A) is fixed in the 
system. In the bulk of the parent phase, the concentration of A (molecules/ 

o 

cm ) is C A (t). At t = 0, it was C A (0). Let us say that at t = 0, the parent 
phase is supersaturated in A so that C A (0) > C A (t). Because of conservation 
of mass, the difference C A (0) - C A (t) is the amount 0 f A deposited in grains 

3 

at time t. The volume of a grain of radius R is (4 tt/ 3)R while the number of 
grains per unit volume with this radius is n(R,t)dR. The average molecular 

* • » • i i i 

volume in a grain is V = X V + XV, where X and X are mole fractions 

A A B B A B 

describing the composition of a grain. The total number of molecules in 

— 1 ' 3 • 

grains of size R is (4 tt/ 3V ) R n(R,t)dR. Of these, X^ are molecules of A. 
Conservation of the total number of molecules of A dictates 


C (0) = C <t) + f dR R 3 n(R,t) (1.4.1) 

A A 3V A o 

where the integral on the right-hand side of Eq. (1.4.1) expresses the number 
of molecules of A in grains of all sizes. 

Since C A (°) > the Parent phase is rejecting A while the second 

phase is concentrating it. The reverse must be true for component B. Hence 

C B (0) < C B (t) ' and 


6 



(1.4.2) 


C B (t) = c b (0) + C~t) x b / dR R3n ( R / fc ) 

3V o 

expresses the conservation of B. 

1 ,5 Grain Size Growth Velocity 

In the previous section, we noted that growing grains of the second 

phase were concentrating A while rejecting B. Hence, in the vicinity of a 

growing grain, there is a net flux of A toward the grain and a net flux of B 

away fran it. Should the grain be dissolving, the situation is just reversed. 

In either case, at every point in the parent phase, it is possible to specify 

local velocities, v A and v R . These velocities are defined in the sense of 

hydrodynamics and express the transport of A and B, respectively, with respect 

$ 

to the laboratory. We may also define a volume averaged velocity, v, given by 

[ 2 ] 


► 


w 


= C V v 
AAA 


C V v 
B B B 


(1.5.1) 


Let n A be the number of molecules of A in the parent phase and let V t be the 

total volume of the system. If is the volume occupied by an A molecule in 

the parent phase, then n is the total volume of A whereas n„V /V = C V is 

AA _ A A t A A 

the volume fraction occupied by A molecules. Likewise C V is the volume 

$ B B 

fraction occupied by B molecules. Hence v is a linear combination of veloci- 
ties v A and v B weighted by respective volume fractions. We may also inter- 
pret by writing on the basis of the discussion above, 

V* ■ 'V» )V S * ‘"bV’b • (1 - 5 - 2 ' 

If in Eq. (1.5.2), the product n A V A were replaced by m A , the mass of A, and 

n B V B replaced by nig, the mass of B, whereas were replaced by the total 

mass, m t = m A + m B' we w °uld readily recognize v* as the velocity of the 

+ 

center of mass. As written, however, it is clear that v in Eq. (1.5*2) is 
the velocity of the "center of volume.” 

Physically, one may indeed see that there must be a velocity associated 
with the center of volume. Suppose < V B and consider a growing grain which 
is concentrating A and rejecting B. A molecules are rushing in toward the 
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grain while B molecules are leaving. Since larger molecules (B) are being 
replaced by smaller molecules (A), the center of volume is shifting; hence, 

v^O. 


It is now useful to define molecular fluxes and by [2] 


*A = °A V A 


(1.5.3) 


H = C v . 

B B B 


(1.5.4) 


± ± 

After adding and subtracting C v and C v from Eqs . (1.5,3) and (1.5.4) , 

A B 

respectively, one obtains 


N = C V* + C (v n - V*) 

A A A A 


(1.5.5) 


N = C v* + C ( v - v + ) . 
B B B B 


(1.5.6) 


By diffusion, we mean motion with respect to some center within the 
parent phase. We choose this reference point to be the center of volume. The 
terms in parentheses on the right-hand sides of Eqs. (1.5.5) and (1.5.6) are 
the diffusive velocities of A and B, respectively. We write the corresponding 
diffusive fluxes 


C R (V ~ V *> = “ D * VC A 

A A A A 


V V B - ''*> ■ - D B • 


Substituting Eqs. (1.5.7) and (1.5.8) into Eqs. 


"a - V* - d a ' c , 


N = Cv - d: vc„ 
B B B B 


(1.5.7) 

(1.5.8) 

(1.5.5) and (1.5.6) gives 

(1.5.9) 
(1.5.10) 


The first berms on the right-hand sides of Eqs. (1.5.9) and (1.5.10) give the 
amount of material transported with respect to laboratory due to the motion of 
the center of volume. This transport is called "advection." The second terms 
on the right-hand sides give the material transported by diffusive motion 


* 


»■ 


■i 


* 
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i + 

measured with respect to the center of volume, where and D 0 are the respec- 
tive diffusion coefficients referred to that center. 

In the remainder of this section, we have two goals. The first is to 

prove that referred to the center of volume, A and B have the same diffusion 

4* -t* 

coefficient, namely D A ~ D B = D , where D is termed the "interdiffusion 
coefficient.” The second is to calculate the rate of diffusion into a growing 
grain, given boundary conditions specifying the concentrations on the surface 
of the grain and in the bulk of the parent phase. To reach these two goals, 
we must first demonstrate some theorems associated with diffusion in this 
frame of reference [2]. 


1.5.1 Two Diffusion Theorems 

lEheorem Is C, VV, + C_ = 0 

A A B B 

Proof : Because the total volume of the parent phase = v t (n A ,n B ) is a 

homogeneous function of first degree, we have by Euler 1 s theorem of thermo- 
dynamics 


V 4 _ = V n + V n 
t A A B B 


(T.1.1) 


The function V^ n A ,n B^ is also a state function (exact differential), so we 
have 



Now form the total differential of Eq. (T.1.1). We find 


(T.1.2) 


■ V"A + V"b + 

Because of Eq. (T.1.2), we conclude from Eq. (T.1.3) that 

n A dV A + = 0 • (T.1.4) 


Since the gradient operator involves just linear combinations of derivatives 
we may write on the basis of Eq. (T.1.4) 



n W + n W = 0 . 
A A B B 


(T.1.5) 


We divide both sides of Eq. (T.1.5) by V t and use C A = n A /V t and = n fi /V t to 


obtain 


C A + C B W B - ° • 


Q.E.D. (T.1.6) 


Theorem 2: V.D* VC. + V_D^ VC_ = 0 

A A A B B B_ _ 

Proof : Multiply Eq. (1.5.7) by and Eq. (1.5.8) by V g and add the results 

-V D* vc - V VC = V C ( v - v*) + V C ( v - v*) 

AAA BB AAA BBB 


= (C V v + C V vj - (VC + VC)f T . 
AAA BBB AA BB 


(T.2.1) 


The first term in parentheses in the second line of Eq. (T.2.1) is by Eq. 
(1.5.1) just v • The second term is 


V A C » + V B C B ' (V a"a + ’ 1 


(T.2.2) 


because 


“ V n + . 

t A A B B 


(T.2.3) 


Using these results, we find for Eq. (T.2.1) 


V D* VC R + V D* VC„ - -v* + v* = 0 
A A A B B B 
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1.5.2 The Interdiffusion Coefficient 




A 


To show that only one diffusion coefficient is required to 
describe interdif fusion, we begin by forming the gradient of both sides of Eq. 
(T.2.2) 


, 0 . „<v s c R * v B c B ) 

- C A '"'a + C B % \ 


where we used the results of Theorem 1. 

Now solve Eq. (1.5.11) for V VC and substitute the result into 

B B 

Eq. (T.2.4). One finds 


0 " ,D » - d b ) v * % 


(1.5.12) 


Equation (1.5.12) can be true in general only if 



where we refer to the common value as the interdiffusion coefficient and 
denote it by D. 


1.5.3 Rate Coefficient and the Effect of Advection 

Consider a system of spherical polar coordinates and a grain of 
radius R with its center located at the origin. The geometry is illustrated 
in Figure 1.1. The fluxes of A and B molecules in the vicinity of this grain 
are given by Eqs. (1.5.9) and (1.5.10). Expressed in spherical polar coordi- 
nates , the first of these is 

3C 

N = C - D — — — . (1.5.14) 

A A 3r 
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* 


* 


Figure 1*1. Geometry for calculating the diffusion into a 
spherical grain of radius , R. 

In writing Eg. (1.5.14)/ we have made note of the fact that Figure 1.1 is 

$ + 

spherically symmetric and that v has only one component/ v / which is in the 
radial direction. The term / C^v , represents radial advection while the 
term, -DSC^/Br/ expresses radial diffusion. Transport in the two angular 
directions is prohibited by the symmetry. We also have used the results of 

A 

Section 1.5.2 and set = D. 

A 

We shall assume that the grain of radius R is isolated from all 
other grains and that there are no sinks nor sources for molecules anywhere 
except on its surface. Thus/ molecules which cross the surface at R must also 
cross the surface of the larger sphere of radius r (broken circle in Figure 
1.1). The total number of molecules of A crossing the sphere of radius r is 
I A and is given by 

3C 

I = 4Trr 2 N = 4Trr 2 C v* - 4Trr 2 D — ^ . (1.5.15) 

AAA 3r 

Because of mass conservation/ I A is a constant independent of r. Now in Eq. 
(1.5.15) replace by its definition given by Eq. (1.5.1). The result is 
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3C 

I = (4irr 2 (C v )V + 4^(0 v )V )C - 4-ni 2 D — - 
A AAA BBBA 3r 


3C 

= + ^JC. - 47rr 2 D — - 

A A B B A 3r 


(1.5.16) 


where we have used Eqs. (1.5.3)/ (1.5.4), and (1.5.15). We have noted that 


I s - 4„r n b 


(1.5.17) 


Our agrument that I A is constant also applies to 1^. This permits us to 
define a constant rate of flow of volume, u, by 


u *= IV + I V 
A A B B 


(1.5.18) 


Using Eq. (1.5.18), we may cast Eq. (1.5.16) in the form of an ordinary dif- 
ferential equation of first order to be satisfied by C A , namely. 




- f ~ u V = - - A - 

3r _ 2 ] A _ 2 # 

47rDr 47rDr 

The general solution to Eq. (1.5.19) is 


(1.5. 19) 


c A ( r ) = (I a A 0 + w exp(-u/47rDr ) 


( 1.5.20) 


where W is a constant of integration. 

By a diffusion-controlled precipitation, we mean one where the 
material just outside a grain is always in equilibrium with the material 
inside. The equilibrium concentration of A just outside the grain is given by 
Eq. (1.2.11). Hence, replacing C^( r) on the left-hand side of Eq. (1.5.20) by 
C a (R) and setting r = R on the right-hand side, we find 


W = exp(u/47iDR) [C (R) - (I / U)] . 


(1.5.21) 


Replacing W in Eq. (1.5.20) by Eq. (1.5.21), we obtain 


C A (r) “ AT + fC A (R) • expf 4^i - 7 ]] * 


(1.5.22) 
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Infinitely far from the grain, the concentration of A is undis- 
turbed by the formation of the grain and has the value, Letting 

r-*oo on the right-hand side of Eq. (1,5.22) and replacing C A (r) by on 

the left, we obtain an equation for I A 

I A = u [C a (R) - c A bUlk ] fexp(u/47rDR) - 1 ] 1 

Define the dimensionless parameter, b, by 

b = u/4ttDR 

then Eq. (1.5.23) may be written 

I A = 4ttDR[C a (R) - C A bulk l b[e b - 1 ] _1 . (1.5.25) 

Since b depends on R, Eq. (1.5.25) is a rather complicated func- 
tion of R to use in further analyses. We note, however, that near the end of 
the precipiation (a circumstance which shall be of interest to us) , the net 
rate of growth given by u diminishes to the point where u/4ttDR << 1. Letting 
e b a 1 + b in Eq. (1.5.25), we find in this limit 

I A = 4ttDR[C a (R) - C A bulk ] . (1.5.26) 

2—1 

Since D has the units of cm sec and R the units of cm, the product, 

3 

4ttDR, has units associated with a bimolecular rate constant, namely, cm 
sec~^. The two species reacting are A and the grains. The rate constant 
being a function of R, we write it as K(R) , where 

K(R) = 4ttDR . (1.5.27) 

An expression for the transport of B, which is the analog of Eq. 
(1.5.26), can be obtained simply by replacing the subscripts A everywhere in 
that equation by B. Note also that when C A (R) < C A bu ** k in Eq. (1.5.26), 
transport of A is in the negative r direction toward the grain. When C A (R) > 
C A bulk , it is in the positive r direction away from the grain. 


(1.5.23) i 


(1.5.24) 
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Because we have made the approximation b << 1 , Eq. (1.5.27) has 
the same form as is found in the absence of advection [7]. The consideration 
of advection r nevertheless , has forced us to pay careful attention to the 
definition of the diffusion coefficients of species A and B. We have found 
that only one diffusion coefficient enters the theory and that it is the one 
governing interdiffusion. In the case that the parent phase is a liquid, D 
and its concentration dependence may be measured directly using the diaphragm 
cell technique [8,9] • 

1.5.4 Diffusion-Controlled Growth Velocity 

The rate of decrease of volume in the parent phase is 
I.V + I V or u as defined by Eq. (1.5.18). The rate of increase of volume 

A A • _ • i * i 

of the second phase grain is I V + IV. If V ± V and V * V then as the 

A A B B AA BB 

precipitation proceeds there is a net change in the overall volume of the 
system. This introduces a new mechanism for growth in addition to diffusion 
and advection. Consider a second phase grain which is expanding faster than 

• t 

the parent phase is contracting, i*e. (I V + I.V.) > (I v + I_ V_) . In this 

AA BB AA BB 

circumstance, the surface of the grain engulfs molecules of the parent phase 
as it advances [10]. This, however, is not a very important effect and shall 
be ignored in the remainder of this report* 

We have been assuming that the second phase grain maintains its 
equilibrium composition during growth. The total rate of arrival of molecules 
to the grain is X A + X B ; hence, for the grain to have its equilibrium composi- 
tion, the molecular rates, I A and I B , must satisfy the equations. 


i 



X A + X B 


(1.5.28) 


I 



*A + *B 


(1.5.29) 


3 

The volume of the second phase grain is (4 tt/ 3)R , and the rate of growth of 
this volume is 


If we multiply and divide the right-hand side of Eq. (1.5.30) through by I A + 
I R and use Eqs. (1.5.28) and (1.5.29), we obtain 


af f*? r3 ) ■ - (i a + v * w 


(1.5.31) 


Using Eq. (1.5.28) again, we may eliminate I 0 from the right-hand side of Eq. 
(1.5.31) and obtain 


dt 


d (il R 3 } - 


-X. 


I • I I 

(x V + XV) 

■ A A B B ■ 


(1.5.32) 


Finally, we compute the indicated differentiation on the left-hand side of Eq. 
(1.5.32) and replace I A on the right-hand side by Eq. (1.5.26), and obtain 


i i i * 


dR (Va + V.) D bulk 

3t ? R ‘ ( C A - C A !R, 1 


Introducing the average molecular volume, V , 


» i i • i 

V “ + X r> V T, 

A A B B 


Eq. (1.5.33) becomes 


(1.5.33) 


(1.5.34) 


dR 

dt 


i 

V D 

V 


( C ; 


bulk 


c a (R>) 


V(R,t) . 


(1.5.35) 


Equation (1.5.35) is the grain size growth velocity, V(R,t) , first postulated 
in Eq. (1.3.1). The next chapter will be devoted to generalizing Eq. (1.5.35) 
to take into account competitive effects due to neighboring grains. 


II. EFFECT OF INTER-GRAIN COMPETITION ON THE GRAIN RADIUS GROWTH VELOCITY 

2.1 Macroscopic Equation of Motion and Its Solution 

We wish to calculate the rate of disappearance of solute molecules 
diffusing freely amongst some randomly placed spherical precipitate grains. 

The grains are stationary, and they and the diffusing molecules are immersed 
in a continuous medium in which the molecules have diffusion, coefficient, 

D. Suppressing the subscripts, A and B, introduced in Chapter I , we let 
C(r,t) be the concentration of diffusing molecules at the position, r, and 
time, t. At the surface of a precipitate grain of radius R, the solute con- 
centration is determined by the Gibbs-Kelvin equation, which takes the form, 

C(R) * C(oo) exp(2yVVkTR) . (2.1.1) 

Whether a particular grain is absorbing solute (growing) or emitting solute 
(dissolving) depends upon the difference (C(r,t) - C(R)). If positive, the 
grain is growing? if negative, it is dissolving. The specific rate of growth 
(dissolution) is K(R) , so that the rate is K(R)(C(r,t) - C(R)). Let P(R)dR be 
the probability that a grain has a radius between R and R + dR. For the 
medium as a whole the rate of change of solute molecules due to the action of 

oo 

the grains is f dRP (R)K(R) (C(r,t) - C(R)). 
o 

Using this integral to take into account the loss of material by incor- 
poration into grains, Fick's second law for C(r,t) is 


v 


3C(r,t) 

at 


oo 

- DV 2 C(r,t) - <|>(r) - / dRP (R)K(R) (C(r,t) - C(R)) 

o 


( 2 . 1 . 2 ) 


where the function <j>(r) has been introduced to represent any continuously 
distributed sources of solute molecules. From a mathematical point of view, 
the presence of <{>(r) permits Eq. (2.1.2) to be treated by Green's function 
techniques. The relevant parts of this and other techniques of mathematical 
physics have been developed in Appendices A through D for those readers who 
may not be familiar with them. It is recommended that those appendices be 
studied now before going further. 
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The physical meaning of <j>(r) depends upon the situation, but one might 
think of <f>(r) as being equal to the quantum yield times the rate of absorbance 
of light producing some photochemical product, which is subsequently precipi- 
tating. Our results will be wholely independent of the form of <J>(r) , however. 

In the steady state 9C(r,t)/3t = 0 and Eq. (2.1.2) becomes 
2 

-DV C(r) = <J) ( r) - / dRP(R)K(R) (C(r) - C(R)) . (2.1.3) 

o 

We may Fourier transform both sides of Eq. (2.1.3) using Eq. (D.26) of 

2 

Appendix D to compute the transform of V C(r) and Eq. (D.10) with k' = 0 to 

00 

compute the transform of \ dRP(R)C(R) which is independent of r. The result 

o 

is 

00 

k 2 DC(k) = 4> ( Jc) - f dRP(R)K(R) (C(k) - (2ir ) 3 C( R) <5 (k) ) . (2.1.4) 

o 


Define 


K T - / dRP ( R ) K ( R ) 


(2.1.5) 


and write Eq. (2.1.4) in the form 


(k 2 D + K T )C(k) = cj>(k) + / dRP(R)K(R)C(R) (2-ir) 3 6(k) 

o 


( 2 . 1 . 6 ) 


where 6(k) is the Dirac delta function discussed in Appendices A and D. 

According to Appendix C and Eq. (D.35) of Appendix D, the Fourier transform of 

2 

the Green’s function G(r) for the operator, DV , is 


G(k) = -1/k 2 D . 


(2.1.7) 


Equation (2.1.7) permits Eq. (2.1.6) to be written 


( 2 . 1 . 8 ) 


00 

/ dRP(R)K(R)C(R)(2u) 3 6(k) ] . 

o 


We next expand the factor fl - G(k)K T ) 


-1 


in a geometric series 


^ = 1 + G(k)K + (G(k)K, l 2 + (G(k)K, J 3 + ... (2.1.9) 

1 - G(k)K T 

and substitute into Eq. (2, 1.8) to obtain 

-C(k) - [1 + (G(k)K T ) + (G(k)K T ) 2 + (G(k)K^) 3 + ...1 

00 

[G(k)J(k) + G(k) / dRP(R)K(R)C(R)(2Tr) 3 6(k)] . (2.1.10) 

o 

We multiply out the brackets in Eq. (2.1.10) collecting together like powers 
of and K(R), since being proportional to the rate of precipitation , they 
are of the same order of magnitude. The result is 

00 

-C(k) = G(k)<J)(k) + [ G ( k) k) <j> ( k) + G(k) f dRP(R)K(R)C(R) (2Tr) 3 6(k) 1 

o 

00 

+ [G(k)K T G(k)K T G(k)^(k) + G(k)K T G(k) f dRP(R)K(R)C(R) (2it) 3 6(k) 1 

o 

+ [G(k)K T G(k)K T G(k)K T G(k)4.(k) + G(k)K T G(k)K T G(k) 

00 

/ dRP(R)K(R)C(R) (2ir) 3 6(k) 1 + ... (2.1.11) 

o 

We next expand K(R) in a series analagous to a cluster expansion in the 
theory of the equation of state of real gases [11]/ 

N 

K(R) = E K( d| R) , (2.1.12) 

d=1 


where N is the number of grains. The term, K( d| R) , represents the contribu- 
tion to K(R) of an ensemble of grains, n d" in number, all competing for a 
diffusing solute molecule. Likewise, Kj, which is the average of Eq. (2.1.12) 
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over the distribution of grain sizes may be developed in a cluster expansion 


as 


N 

K = l K (d) (2.1,13) 

d=1 

where 

00 

K T (d) = / dRP(R)K(d| R) . (2.1.14) 

o 

We put Eqs. (2.1.12) and (2.1.13) into Eg. (2.1.11). For simplicity of nota- 
tion, we suppress everywhere the dependence upon the Fourier transform vari- 
able, k, and write 


A A A A AAA 

-C = G(J) + [G(K t (1) + K t (2) + K t (3) + ...) G(J> + G 

°° 3 

f dRP(R) (K{ 1 | R) + K( 2 1 R) + K(3 |r) + ...) C(R)(2tt) 6] 

o 

+ [G(K t (1) + K t (2) + K t (3) + ...) G(K t (1) + K t (2) + K t (3) + ...) 

gJ + G(K t (1) + K t (2) + K t (3) + ...) G 
00 

/ dRP(R) (K(l|R) + K(2 |r) + K(3|r) + ...) C(R)(2t;) 3 6] (2.1.15) 

o 

+ [G(K t (1) + K t (2) + K t (3) + ...) G(K t (1) + K t (2) + K t (3) + ..,) 

G(K t (1) + K t (2) + K t (3) + ...) gJ 

+ G(K t (1) + K t (2) + K t (3) + ...) G(K t (1) + K t (2) + ^(3) + ...) 

oo 

G / dRP(R) (K( 1 | R) + K( 2 1 R) + K(3 |r) + ...) C(R)(2tt) 3 6] . 
o 

On the right-hand side of Eq. (2.1.15), terms of the same order need to be 
collected together. Note, for example, that Krp(2) and (K,p(1) are of the same 
(second) order, while (K^l)) 3 , ^(1)^(2), and *^(3) are of the same (third) 
order. The result is 
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-C(k) = G(k) <J>(k) + [G(k)K T (1)G(k) J(k)] + [G(k)K T (2)G(k)J(k) 

+ G(k)K T (1)G(k)K T (1)G(k)J(k)] 

+ [G(k)K T (3)G(k) J(k) + G(k)K T (1)G(k)K T (2)G(k)^(k) 

+ G(k)K T (2)G(k)K T (1)G(k)<f,(k) + G(k)K T (1)G(k)K T (1)G(k)K T (1)G(k) J(k)] 


+ • • • 


( 2 . 1 . 16 ) 


+ [G(k) / dRP(R)K(3|R)C(R) ( 2 TT ) 3 6 ( k) + G(k) f dR*P(R' )K(1 |r' )G(k) / 


dRP(R)K(2|R)C(R)(2ir) 3 5(k) 


+ G(k) / dR , P(R')K(2|R')G(k) J dRP(R)K( 1 | R)C(R) (2u) 3 6(k) 

o o 

A “ A 00 A °° 

+ G(k) / dR"P ( R" ) K ( 1 | R" ) G ( k) f dR'P(R' )K( 1 |r' )G( k) / 


dRP(R)K(l|R)c(R)(2ir) 3 6(k)] 


+ 


• • • 


This completes the analysis based upon a macroscopic point of view. In 
the next section (Section 2.2), we reanalyze this problem, working on the 
assumption that the locations of all the grains are known. Subsequently, in 
Section 2.3, we shall average the result from Section 2.2 over a distribution 
of grain sizes and locations (configurations) to obtain an equation for com- 
parison with Eq. (2.1.16). This comparison will be carried out in Section 2.3 
and will result in a revision of the formula, K(R) = 4 ttDR, to take into 
account the effect of grain density on the diffusion process. 

2.2 Microscopic Equation of Motion and Its Solution 

In this section, we seek to represent on a microscopic scale the diffu- 
sion of molecules among the grains. We specifically consider a continuous 
parent phase in which some grains of second phase are imbeded. The grains are 
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N in number and have centers located by the set of coordinates = (r.j,r ? , 

..., r^) • let c(r?{r^}) be the concentration of molecules diffusing among the 
particular configuration , {r^}, of the grains* The equation of motion for 
c(r;{r^}) is Fields second law, which in steady state takes the form 

2 N 

DV c(r) = - E q.6(r - r. ) (2,2.1) 

i=1 1 1 

where we have intentionally suppressed in c(r) the dependence upon {r^}. The 
i-th grain is located at r^ , has radius and "strength” q^, which expresses 

its ability to absorb or emit molecules. The spatial location of this 
strength is given by 6(r - r^) • Inasmuch as <5(r - r^) represents the location 
of a point particle, it is in seeming contradiction with our assignment of a 
radius, , to the grain. This concept will be dealt with carefully when we 
determine the values of the {q^}. For illustration, the geometry of just two 
grains is shown in Figure 2.1. More on the Dirac delta function can be found 
in Appendix A, whereas Fick’s laws of diffusion are discussed in Appendix B. 


Z 
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Figure 2.1. Coordinates locating the centers of two grains of radii 

R 1 and R 2 * respectively. 


Each of the N grains is assumed to be in chemical equilibrium with the 
parent phase just outside its boundary. Hence, the concentration just outside 
each grain is specified by the Gibbs-Kelvin equation. Every point on the 
surface of a grain, in fact, must satisfy the Gibbs-Kelvin equation, so we 
have 

C(r i + = C(R i } = C(oo) ex p( 2yV/kTR i ) , i = 1,2, ..., M (2.2.2) 

where - l R il- If the grains are sufficiently dilute, then R^ = I I « 
|r^|, and we may to sufficient accuracy replace Eq. (2.2.2) by 

c ( r i ; { r i } ) - C(R ± ) » C(-) exp(2yV/kTR i ), i = 1,2, ... M (2.2.3) 

Eqs . (2.2.3) may be used to determine the {q^} as follows. 

2 

As shown in Appendix C, the Green's function for the operator DV is 
given by 

G(rjr') (2.2.4) 

47rD|r - r' | 


which satisfies the boundary conditions 
lim G(r|r' ) = lim VG(r|r' ) = 0 . 

r>oo r-^oo 


(2.2.5) 


Hence, according to Appendix C, we may apply Green's theorem in the form 

c(r) = f d 3 r' G(r|r' )f(r' ) . (2.2.6) 

If on the right-hand side of Eq. (2.2.1) we replace i by j and make the 
identification 


N 

f(r) = - 4> (r) - r q .6 ( r - r.) (2.2.7) 

3-1 D 3 
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then we may solve Eq. (2,2. 1) for c(r) in the form 


N 


-c 


(r) = f d r G(r|r' ) [<f>{r* ) + E q.6(r' - r.)] 

j-1 3 3 


N 


» / dr G(r| r' ) <t>(r' ) + E q.G(r|r ) 

j=i 3 3 


( 2 . 2 . 8 ) 


where we have used Eq. (A. 9) of Appendix A to evaluate the integral involving 
the delta function. 

We next use Eq. (2.2.8) to evaluate c(r) on the boundary of the i-th 
grain. We note that this value, C(R^) is established hy Eq. (2.2.3). We find 


-CdO 


- / d r'Gfi^ 


r.) 

D 


(2.2.9) 


*i 


To obtain Eq. (2.2.9) from Eq. (2.2.8), care has been exercised in 
evaluating G(r|r_.) on the surface of each sphere. The cases i=j and i*j must 
be distinguished. To evaluate G(r|r_.) on the surface of the i-th sphere where 
i=j, we have set r = r^ + R, , and written 


G(r. + R. | r . ) - = G(R ). (2.2.10) 

1 11 4irD|r + \ ~ r J 4 ttD | R^ | 4 ir DR i 


This accounts for the term q^Gd^) on the right-hand 
obtain G(r|r\) on the surface of the i-th sphere (i* 
, and written 


side of Eq. (2.2.9). 
j) , we have set r = 


To 


G(r. + R. |r .) a = G(r. - r.) (2.2.11) 

1 1 5 4.d I r ± . - r 3 | 4.D | r ± - r.| 1 3 


Eq. (2.2.11) follows because the dilute nature of the sinks which implies 
| R ± | << |r^ - r^| and |r^ + R^ - r^| = jr^ - r j|* Eq. (2.2.11) accounts for 
all the terms within the sum on the right-hand side of Eq. (2.2.9). 
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By virtue of Appendix D and the fact that G(r, |r' ) = -(4 ttD)~ 

|r^ - r 1 | = G(r - r|), we may recognize the integral term on the right-hand 

side of Eq. (2. 2.9) as a convolution and write 


G*<J>(r) = / d 3 r'G(r - r* )<f>(r') . 


( 2 . 2 . 12 ) 


Substituting Eqs. (2.2.10) and (2.2.12), we solve Eq. (2,2. 9) for q i , as 

N 

q = 4 ttDR [G*(j,(r ) + z G(r. .) q. + C(R )] (2.2.13) 

ii <j D 3 i 

*i 


where we write G(r. - r .) = G(r, .) = -(47rDr. .) -1 with r . . = I r , .1 = 

i 3 ij ij 1 3 ±3' 



As it stands, Eq. (2.2.13) does not determine q^ because the equation 

depends for its evaluation upon a knowledge of the {q.}, j*i. By the method 

1 3 

of successive approximations, however, we may use Eq. (2.2.13) to develop Eq. 
(2.2.8) in a power series. To do this, we replace j by i in Eq. (2.2.8) and 
write 


N 


-c(r) = G*(f>(r) + E 
i=1 




(2.2.14) 


Now substitute Eq. (2.2.13) into the right-hand side of Eq. (2.2.14). The 
result is 

N 

-c(r) = G*4>(r) + E G(r - r ) 4 ttDR [G**(r ) + 
i=1 ill 

N 

E G(r ) q s + C( R. ) ] . (2.2.15) 

j=1 3 3 


We expand Eq. (2.2.15) further. In Eq. (2.2.13), replace i by j and j by £ to 
obtain 
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(2.2. 16) 


q 


j 


4ttDR 


+ t i G< v q * + c<R j > ' • 


Substitute Eq. (2*2.16) for qj on the right-hand side of Eq. (2*2.15) to 
obtain 


N 

-c(r) = G*(J,(r) + E G(r - r^TrDR fG^frp + CfR^l 
i=1 

N N 

+ E £ G(r - r. )4TrDR.G(r. ,)4 ttDR .[G*<(,(r .) + C(R.)1 (2.2.17) 

i-1 j-1 1 1 13 31 3 3 

*i 


N 

+ E 
i=1 


N N 

E E G(r - r, )4 ttDR G(r. . )4 ttDR .G( r ) 4 ttDR Tg*^ ( r ) +C(R )} . 
. 1 i 13 3 3l £ L Y l l 


In Eq. (2.2.17) in the triple sum on the right, we have replaced q^ by 
4 ttDR^[g*<() ( r^) + C(R^)] which is the leading turn in the next substition. This 
"closure" step is a necessary approximation if we are to have an equation for 
c(r) which is free of explicit dependence upon any of the fq^}. Equation 

(2.2.17) is the result of only two successive approximations? the procedure, 
however, may be carried out to any higher degree of approximation desired. 

It is now worthwhile to examine a product such as G(r - r^) 4 ttDR^ 

G*<j>(r^) which appears, for example, as a factor within the single sum in Eq. 

(2.2.17) . If we define the "scattering" operator T i (r,r* ) by 


(r,r* ) = 4 ttDR^ 6 (r - r i )6(r l - r ) 
then the desired product can be written 


(2.2.18) 


G(r - ) tt 4DR^G*(J) ( r ± ) 

= G(r - r 1 )4 7 rDR i / cPr'Gfi^- rp<f,(r’ ) (2.2.19) 

= / d 3 r* "G(r - r' ") f d 3 r' '4irDR 6(r' • ' - r )S(r" - x^) 

/ d 3 r , G(r' ' - r* ) 4> (r* ) = G*T ± *G*<), ( r) . 
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Using Eqs. (2.2.18) and (2.2.19), we may write Eq. (2.2.17) in the form 

N 

-c(r) = G*,j>(r) + E G*T *[G*$(r) +C(R.)] 
i =1 1 


N N 

+ E E G*T. *G*T .*[G*(f)(r) +C(R.)] 

i =1 j =1 1 3 ° 

* i 

N N N 

+ E E E G*T. *G*T ,*G*T *[G*<j>(r) + C(R )] + ... (2.2.20) 

i =1 j =1 1=1 1 3 4 

*j 

The ultimate goal of this section is to obtain a result which can be 
compared with Eq. (2.1.16) of the previous section. In order to do this, we 
form the Fourier transform of both sides of Eq. (2.2.20). The transform of 

a 

-c(r) is -c(k), but the transform of the right-hand side must be analyzed with 
some care. Because of the use of the convolution symbol throughout Eq. 
( 2 . 2 . 20 ), it would seen that the transform of each term within the various 
sums would simply be the product of the individual transforms. However, Eq. 
(2.2.18) involves the product of delta functions and is not explicitly a 
function of r - r' , which was a property exploited in going from Eq. (D.14) to 
(D. 15) of Appendix D. We will see, nevertheless, that Eqs. (D.18) and (D. 19) 
apply if we define carefully the transform of T^r/r' ). Consider the trans- 
form of G*T i *G*<(>(r) . Since this is a transform of a convolution, we denote 
it by 


(FTC) = / d 3 r e _lk * r G*T ± *G*<J>( r) 

= Jd 3 r e" ik * r /d 3 r‘ • ’G(r-r" ') f d 3 r' 'T^r' • • ,r’ • ) fd 3 r'G(r' '-r' ) 4 , ( r * ) 

= /d 3 r/d 3 r' "/d 3 r" fd 3 r' e“ ik ’ r G( r-r* " )T (r" * ,r* • ) G(r* ’-r' ) «j>(r' ) . (2.2.21) 

3 3 

Replace the variable r by £ = r - r" ' , where d r = d 5 , so that Eq. (2.2.21) 
reads 
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(FTC) = / d 3 ge‘ ik * ? G(£) / d 3 r* 1 1 

/ d 3 r d 3 r , e" ik,r ' "^(r' • ’ ,r* *)G(r* '-r* ) + (r' ) . (2.2.22) 


Replace r' ' ' by p = r' ' '-r* ' , where d 3 r’ ' ' = d 3 p, and for T^(r' ' ' ,r* ') 
substitute 


(r* 1 1 ,r* 1 ) = 4^^DR i S(r ,,, - r ± )6(T" - r ± ) (2,2.23) 

so that Eq. (2.2.22) reads 

(FTC) = /d 3 £e~ llc *^G( £)4irDR i Jd 3 r' 'e -lk>r ' ' 6(r* * - r )/d 3 pe” ik * p 6(|»-(r - r' • ) ) 

x /d 3 r’G(r' ' - r' ) <t>(r* ) 

= /d 3 £e~ lk *^G( g)4irDR i j’d 3 r' , e -lk>r ' ' 6(r" - r ± ) e’ ik ' (r i* r ’ ’ ) 

x /d 3 r’G(r" - r’ ^(r* ) . (2.2.24) 


In Eq. (2.2.24) , we have used the right-hand side of Eq. (D.7) to evaluate the 
integral over d 3 p. 

For 5(r* * - r^), it is convenient to use the Fourier representation 
given by Eq. (D.8) which reads 

6 (r' • - r t ) = 6 ( 1 ^ - r' • ) = (1/2n ) 3 / d 3 k' e lk ' * (r i r ' ' ) . (2.2.25) 

In writing Eq. (2.2.25), we have exploited the fact that the delta function is 
an even function of its argument [see Eq. (A. 14) of Appendix A]. Substitut- 
ing Eq. (2.2.25) into Eq. (2.2.24) yields 
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(FTC) = /d 3 Ce‘ lk ‘ 5 G(5)47tDR i /d 3 r"e" i,c ’ r " (-L] 3 fd 3 k‘ e ik ' * {r i " 
x /d 3 r'G(r' ' - r ' ) 4> ( r* ) 

= [/d 3 Ce‘ ik ’ 5 G(e)][4 1 rDR i (^) 3 / d 3 k'e i(k ' " k) * r i] 
x [/d 3 r' 'e” ik ' * r ' ' J d 3 r'G(r' ' - r' )«j>(r' )] 

= [/d 3 Ce _lk * ? G(5)][4TTDR i (^) 3 /d 3 k , e l(k, " k) * r i P(k|k')l 
x [/ d 3 r' 'e * k * r j d 3 r'G(r* ' - r' )<J>(r' )] . 


r' • 


-ik. (r - r") 

e 1 


(2.2.26) 


(2.2.27) 


Note that within the last bracket on the right-hand side of Eq. (2.2.27) , we 
have replaced k’ by k, and in the next to last bracket introduced the symbol , 
P(k|k' )• The effect of P(k|k' ) is to change k to k' in all functions on its 
right. This permits us to go from Eq. (2.2.26) to (2.2.27) and to identify 

a a 

the last bracket as G(k)<j>(k) using Eqs. (D.18) and (D.19). We further define 

T\(k) = 4irDR i ) 3 /d 3 k' e l(k ' " k) * r i P(k|k’ ) . (2.2.28) 

With Eq. (2.2.28) and the results of this paragraph, we may write Eq. (2.2.27) 
in the form 


(FTC) = GfkjT^kJGfkJ^k) . (2.2.29) 

Note carefully that because of the definition of P^k* ), T. (k) is an operator 

A A 

in k space, and it and functions of k, such as G(k) and <{>(k) , do not commute. 

To complete the calculation of the transform of Eq. (2.2.20), we need 
also to compute the transform of terms of the type G*T^*C(R^). Let the trans- 
form of a term of this type be called FTD. The manipulations for computing 
FTD are the same as Eqs. (2.2.21)-(2#2.27); namely. 


(FTD) = [/d 3 5e _ik ’ ? G(5) 1 |4nDR i f2^) 3 


fdW (k - k,) * r i 


P(k|k* ) ] 


X 



-ik.r 1 1 

e 


C(R. ) 1 . 
1 


(2.2.30) 
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Equations (2. 2. 27) and (2*2.30) are identical except for the last bracket in 
Eq. (2*2.30) involving the constant C(R^) , which replaces G*<f>(r) in Eq. 
(2.2.27). Since C(R^) does not involve r M , it may be factored out of the 
bracket. After factoring, if we then compare the last integral in Eq. 

(2.2.30) with Eq. (D.10) with k* = 0, we have 

C(R i )/a 3 r”e" llc * r "= C(R i )(2 1 r) 3 5(k) . (2.2.31) 

Combining Eqs. (2.2.30) and (2.2.31) and using the previous results of this 
section, we find 

( FTD) = G(k)T i (k)C(R i )(2 1 r) 3 6(lc) . (2.2.32) 

Applying Eqs. (2.2.29) and (2.2.32) to Eq. (2.2.20), we get 
N 

-c(k) = G(k)J(k) + E G(k)T. (k) CG(k) J(k) + C(R. ) (2Tr) 3 6(k) ] 

i=1 1 1 

nn aa 

+ E E G(k)T (k)G(k)T.(k)[G(k)i(k) + C( R . ) ( 2 tt) 3 $ ( k) ] (2.2.33) 

i=1 j=1 3 3 

N N N A A 

+ E E E G ( k) T , ( k) G ( k) T . ( k) G ( k) T ( k) [ G ( k) ^ ( k) + C(R ) (2Tr) 3 S(k) 1 . 

i-1 j=1 £-1 1 3 1 ' 1 

i-*j 3*1 

Equation (2.2.33) applies to a configuration of grains with specified radii 
and fixed locations. Ordinarily, there are so many grains that it is 
impossible to specify either their radii or their locations. Rather, only the 
probability of a given radius and probability of a given location can be 
specified. In Section 2.1, we let P(R)dR be the probability of a grain having 
a radius lying between R and R 4* dR. If the volume of the system is V, the 
probability of finding the center of the i-th grain within an element of 

3 

volume is d r^/V • We now specify that the grains are randomly located and 
exert no long range forces on one another. Furthermore, we shall take their 
concentration small enough so that the fact that they must not overlap may be 
ignored. For these reasons, their positions are uncorrelated and the 
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probability of a given configuration (specified set of values for the {r }) is 

N 3 

just the product of individual probabilities, namely, II (d r./V). In Eq. 

i-1 1 

(2.2.33), c(k) is an implicit function of the radii { } and the positions 
{r.}. After drawing the {R. } from the distribution P(R^), the 
configurational average of c(k) should be equal to C(k) in Eq. (2.1.16). We 
define 

N oo N 

«c(k)>> = n / dR, P(R. ) v" N n / d 3 r. c(k) = C(k) (2.2.34) 

i=1 o 1 1 j=1 3 

A 

where <<c(k)>> symbolizes the average of c(k) over both the distribution of 
grain radii and the distribution of grain configurations. 

We form the average of both sides of Eq. (2.2.33), taking note that 

A 

because of (2.2.28), T^(k) depends upon both R^ and r^. The result is 


A A 


N 


A A 


«c(k)>> = G(k) <)>(k) + E [G(k)<<T. (k)>>G(k) 4>(k) 

i=1 1 


+ G(k) «T i (k)C(R i )»(2 ir ) 3 6(k) ] 


N N 


A A 


+ E E [G(k)<<T. (k)G(k)T . (k) >>G(k) <j>(k) + 
i=1 j=1 1 3 


*i 


G(k)<<T i (k)G(k)T j (k)C(R j )»(2Tr) 3 6(k)l 


N N N 

+ E E 


A A 


AAA A A 


i=1 j=1 i 


E [G(k)<<T. (k)G(k)T ,(k)G(k)T (k)>>G(k) <(.(k) + 
_ ^ D x 


A A A 


A A 


+ G(k)<<T. (k)G(k)T .(k)T .(k)G(k)T (k)C(R )»(2 tt) 6(k) + ...] . (2.2.35) 

1 1 3 • - 
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Note that P(R^) is normalized to unity, namely, 

CO 

f dR P(R ) = 1, £ = 1, 2, ..., N 

o * * 

and also that 

v -1 / d 3 r^ =1, £ = 1, 2, ..., N . 

Hence, averages of typical terms shown on the right-hand side of Eq. 
follow. For example. 


N 

«G( k) <J> ( k) >> = n / 
i=1 o 


dR i P(R i ) 


-N 


N 

n 

j=i 


/d 3 r^. G( k) <j> ( k) = G(k)(()(k) 


A 

because neither G(k) 
typical average is 


A 

nor <jj(k) 


depends upon any of the 



<<T . (k)G(k)T . (k) >> 
i J 


N °o N 

n f dR p(r )v n 

£=1 o 1 1 m=1 


f d 3 r T, (k)G(k)T , (k) 
J m 1 J 




fd 3 r . T. (k)G(k)T . (k) 
J 3 i J 


N 

n 

£=1 

*i 


/ 

o 


d V (R £ )V 


-(N-2) 


N 

n 

m=1 

*i 





(2.2.36) 

(2.2.37) 
(2.2.35) 

(2.2.38) 
Another 


(2.2.39) 


32 


Note that because of Eqs. (2.2.36) and (2.2.37), the integrals involving 
neither i nor j are unity. 

This completes the analysis based upon the microscopic point of view. 


2.3 Rate Coefficient and Diffusion-Controlled Growth Velocity 

By virtue of the connection between <<c(k)>> and C(k) expressed by Eq. 
(2.2.34), Eqs. (2.1.16) and (2.2.35) express the same thing. Both Eqs. 
(2.1.16) and (2.2.35) have been organized so that terms of like order have 
been grouped together. Now comparing terms of like order between the two 
equations, we obtain 

N 

K t (1) - E <<T ± (k) >> (2.3.1) 

i=1 


* N N 

K (2) + K ( 1 )GK ( 1 ) = E E «T. (k)G(k)T . (k)» 
T T T . , . . i i 

1=1 j=1 

* i 


K t (3) + K t (1)GK t (2) + GK t (2)GK t ( 1) + K^dJGK^d) - 


MW A A A A /V 

Z Z « T. (k)G(k)T ,(k)G(k)T (k)>>, 
j=1 1 = 1 1 3 * 


etc. The terras involving the grain surface concentration yield, 

00 N A 

/ dRP(R)K( 1 | R)C( R) = Z «T. (k)C(R )>> 

O i=1 1 1 


/ dRP(R)K( 2 | R) C(R) + f dR'P(R’ )K( 1 |r’ )G / dRP(R)K( 1 |R)C(R) (2.3.4) 

o o o 


N 

Z 

i=1 


(2.3.2) 


(2.3.3) 
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N N A A A 

= E E «T. (k)G(k)T . (k)C(R .)>> (2.3.5) 

i=1 j=1 3 3 

*i 


oo oo oo 

/ dRP (R)K(3 |r)C(R) + / dR'P(R' )K( 1 |r' )G / dRP ( R)K( 2 | R) C(R) 
o o o 


oo oo 

+ / dR'P ( R' )K( 2 | R' )G / dRP(R)K( 1 | R)C(R) 
o o 


oo A oo a oo 

+ / dR"P ( R" ) K ( 1 | R" ) G f dR'P(R' )K( 1 | R' )G / dRP (R)K( 1 | R)C(R) 
o o o 


N N N 

= Z E T, <<T. (k)G(k)T .(k)G(k)T (k)C(R )>> . (2.3.6) 

i=1 j=1 1 111 


Equations ( 2 . 3 • 4)- ( 2 • 3 .6) contain no information not already contained in Eqs. 
(2.3.1 )— (2.3.3) j hence, we shall choose Eqs. (2.3.1 )— (2.3.3) for further 


analysis and proceed to the computation of the averages <T. >, 

and <<T.GT .CT >>. 
i 3 i 


<<T . GT .>> 
1 3 


2.3.1 Computation of <<T. (k)>> and K t (1) 


According to Eq. (2.2.28) and the definition of an ensemble 
average expressed by Eq. (2.2.34), the meaning of <<T^(k)>> is 


«T ± (k)» = 


dR P ( R . )4irDR.V _1 fd 3 r. 
i l l 1 i 

o 

We interchange the order of integration with respect to d 3 ^ and d 3 k' on the 
right-hand side of Eq. (2.3.7) and use Eq. (D.10) to recognize a delta func- 
tion. Because there is no function (i.e., only unity) to the right of 
P(k|k' ) on which it can act, we may set P(k|k' ) = 1 . We find. 






P(k k* ) 


(2.3.7) 
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<<T.(k)>> = \ dR . P ( R . ) 4 ttDR , V f— 1 fd k 1 fd r.e 1 

i J 11 l 1 2 tt ' J J l 

o 

CO 

= f dR. P(R. )4TrDR.v" 1 f-ll 3 fd 3 k' (2Tr) 3 6(k’-k) 
l i i ( 2n 


= 4irDV 1 f dR. P(R. )R. 

J ill 

o 

= 4ttDV _1 <R> 


where we have used Eq. (A# 9) with r — k 1 , r f = k, and f(r) = 1, and 
<R > , we mean 


<R> = f dR. P(R. ) R. 

J l l l 

o 


According to Eqs. (2.3.1) and (2.3.8) 


N N -1 

K (1) = £ «T . (k) >> = £ 4 ttDV <R> = 4ttD<R>(N/V) 

T i . - 

1=1 1=1 


because there are N terms in Eq. (2.3.10), each equal to 4 ttD<R>V \ 


K T (1) = 4ttD p <R> 


where 


p = N/V 


is the number density of grains. 

A A A 

2.3.2 Computation of <<T, (k)G(k)T (k) >> and K T (2) 
According to Eqs. (2.2.28) and (2.2.39) 


«T . (k)G(k)T . (k)>> 
i J 


00 


oo 


= / dR i P(R i ) / 
o o 


dR.P(R.) V _1 fd 3 r.4irDR. (-rl) 3 f 

"1 1 *1 l v 7*ir' - 


(2.3,8) 
where by 


(2.3.9) 


(2.3.10) 


Hence/ 

(2.3.11) 


(2.3.12) 
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i(k' - k).r 

d k'e P(k|k’)G(k) 


xV 1 fd 3 r ,4irDR . f— ^-l 3 fd 3 k’ 'e 
D J 


i(k' ' - k) .r. 


P(k|k") (2.3.13) 


- 2 , 


= / dR i P(R i ) / dR _.P ( R ^ ) V (4irDR i ) (4ttDR ^) 


o * . •, i(k - k’ ).r 

x / d k'G(k)(— ) / d 


, . o , i( k' ' - k' ).r. 

x /d 3 k»(^l) 3 /d 3 r.e 3 


where we have used the definition of P(k' '(k 1 ) to exchange k for 

i(k' f - k).r 

k' in G(k) and in e 

We have also interchanged the order of integration with respect 
to d 3 r^ and d 3 k* and with respect to d 3 r. and d 3 k" . In Eq. (2.3.13), the 
integral over d r^ is equal to (2 tt) dfk’-k) while the integral over d r . is 

3 

equal to (2 tt) ). After these identifications, we note that 

/d 3 k'G(k' )6(k'-k) = G(k) by Eq. (D.12) and /d 3 k' , 6(k"-k , ) = 1 by Eq. (D.11). 
Equation (2.3.13) becomes 


<<T i (k)G(k)T_. (k)>> 


= V _2 (4ttD) / 
o 


dR i P(R i )R i (4 1T D) / 
o 


dR..P(R .)R,G(k) 


= v" 2 (4ttD<R>) (47rD<R>)G(k) 


(2.3.14) 
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We observe that on the right-hand side of Eq. (2*3.2) there is a 
double siim, the inner sum having one less term (j*i) than the outer. The 
outer sum (over i) has N terms (i= 1, • • • , N) while the inner sum (over j) 
has (N-1) terms (j = 1, • i - 1, i + 1, •••/ N). Hence, the double sum has 

A. 

N (N— 1 ) terms. Because G(k) on the right-hand side of Eq. (2.3.14) depends 
upon neither i nor j, each of these terms is identical, and Eq. (2.3.2) 
becomes 


K t ( 2) + K T (1)G(k)K T (1) = f^l(~l(4 1 TD<R>) 2 G(k) 

= (4ttD p <R>) 2 G(k) . (2.3.15) 

To get Eq. (2.3.15), we have taken the "thermodynamic limit," wherein 
N-m», V-m» while the ratio, p = N/V, is held constant. When we substitute Eq. 
(2.3.11) into the left-hand side of Eq. (2.3.15), we find 


K,p( 2 ) = 0. (2.3.16) 

2.3.3 Computation of «T ± (k)G(k)T . (k)G(1c)T^ (k)>> , K T (3), and 
Higher Order Terms. Use of Diagrams 

A A A A A 

In computing <<T. (k)G(k)T . (k)G(k)T (k)>> , we need to distinguish 
i j 

two cases. 

Case Is i*£ 

In this case, we have 

<<T . (k)G(k)T ,(k)G(k)T (k)>> 
i 3 Jc 



dR i P(R i ) / 
o 



3r i 4 * D V^ 3 r 


i(k'-k).r 

d k’e P(k|k' )G(k) 


x V 1 /d 3 r j 4irDR j ( I i) 3 


i(k"-k).r. 
/ dx* * e 3 


P(k| k' * )G(k) 


it ITT i(k'"-k).r 

x V" 1 /d 3 r M * ,(~) 3 / d 3 k' ' ' e 


P(k k' ' ' ) 


(2.3.17) 
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To simplify Eq. (2.3.17), we make repeated use of the techniques employed to 
go fran Eq. (2.3.13) to Eq. (2.3.15). The result is 


<<T. (k)G(k)T ,(k)G(k)T (k)>> = V -3 (4 ttD<R>) 3 (G( k) ) 2 
1 3 ** 


(2.3.18) 


Case II: i=£ 

Evaluation of the triple sum on the right-hand side of Eq. 

(2.3.3) requires j*i and l* j. Not ruled out, however, by these restrictions 
is the possibility, £=i. We now treat this case. Because £=i, we must return 
to Eq. (2.2.39) to reexamine the averaging process. On the basis of Eq. 
(2.2.38), an average over three |T_,(k)} operators is 

AA/tAA. N 00 N 

«T. (k)G(k)T ,(k)G(k)T. (k)>> = II / dR P(R)V~ N II f 

1 j 1 „ . 1 i a . - 

1= 1 o m=1 

d 3 r T. (k)G(k)T , (k)G(k)T. (k) 
mi j 1 



/d 3 r.T (k)G(k)T.(k)G(k)T (k) n 
3 3 £=1 


*i 


f 

o 


dR P(R ) 
i i 


N 

n 

m=1 

ti 



r 

m 



/ 


d 3 r^T’ i (k)G(k)T^ (k)G(k)T ± (k) 


(2.3.19) 


where we used Eqs. (2.2.36) and (2.2.37) to evaluate integrals not involving i 
and j. Substituting Eq. (2.2.28) into Eq. (2.3.19) we obtain 



<<T i (k)G(k)T j (k)G(k)T i (k)>> 



dR.jPOO / dR]P(R|) / dR P(R ) V’ 1 / d 3 r i 4 1T DR i (-^ : ) 3 / 
o o 


d 3 k' e 


i(k'-k) .r. 


1 P(k|k’)G(k) 


x V 


-1 


3 

d r . 4 ttDR , 
3 3 


(aiV 


f d 3 k' • 


i(k"-k).r 

e J P(k|k" )G(k) 


x 4 ttDR^ 



i(k* ' '-k) .r ± 
e 


p(k|k’ ••) . 


(2.3.20) 


In Eq. (2*3.20) , we gather together the integrals over dR^, 
dR^ / dR_. , compute them, and obtain the result (4 ttD<R>) • Operating on unity, 
P(k|k* '*) has no effect. The remaining interchange operators must be care- 
fully computed, however. We write the relevant factors fr om Eq. (2.3.20), 

i(k' *-k) «r i(k* 1 *-k) *r 

P(k|k')G(k)e 3p(k|k*)S(k)e 1 P(k|k"’> 


A i(k' '-k).r, A i(k' "-k")T 

= P(k|k' )G(k)e 3 G(k")e 1 


i(k"-k')T. A i(k" , -k’ ').r 

= G(k')e 3 G(k")e 1 . (2.3.21) 


Taking note of Eq. (2.3.21) and the other results of this paragraph, we write 
Eq. (2.3.20) in the form 
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A A 


A A 


<<T . (k)G(k)T . (k)G(k)T. (k) >> 
i 3 i 


I « I 


= (4ttD) / dR i P(R i )R i (4i T D) / dR i P(R i )R i (4 1 rD) f dR j P(R j )R_.v" 2 f^) 3 f 


d 3 k’G(k' ) f— ^-) 3 fd 3 k' 'G(k' ' ) f —^') 3 fd 3 k’ • ’ 
z-n J ' 


2tt 


r , ,3 -ir. .(k-k’+k 1 , -k , ") nrf .3 '^.(k’-k") 

* [/d V i ][/dr« ] 


(2.3.22) 


In Eq. (2.3.22) 
( 2tt ) 3 6 (k-k’+k* • 


, the first integral in brackets is equal to 

-k"’) = (2tt) 3 6(k» "-k'+k"-k) = (2tt) 3 6 (k* "-(k'-k' '+k) ) 


because the delta function is even* The second integral in brackets is 
(27T) 3 5(k , -k M ) • Putting these results in Eq* (2*3.22) , we obtain 


«T i (k)G(k)T j (k)G(k)T i (k)» 


= ( 4ttD<R> ) 3 v” 2 p2~) 3 f d 3 k'G(k') 


X 


/d 3 k' ' 6(k' '-k' )G(k' ' )/d 3 k' ' ' 


6(k' ' '-(k’-k' '+k) ) 


= (4ttD<R>) 3 v” 2 ( 2 “) 3 / d 3 k'(G(k')) 2 . 


(2.3.23) 


In obtaining the last line of Eq. (2.3*23), we have noted that the integral 

3 3 

over d k''' is unity due to Eq. (D.11) while the integral over d k' 1 is 

G(k* ) because of Eq. (D.12). To complete the evaluation of Eq. (2.3.23), we 

need to compute the remaining integral over d 3 k* • We may replace 

a. O -a 

G(k* ) by ( - 1 /D ( k * ) ) using Eq. (D.35). Although written G(k* ) , this function 

in fact depends only upon the length of the vector, k* [see Eq. (D.35)] , and 

not its direction. Hence, it is worthwhile to compute the integral over d 3 k* 

using polar coordinates (k 1 ,0 , <j> , ) as follows 

K , K 
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«T i (k)G(k)T j (k)G(k)T i (k)» 

oo 2 tt 7T 

= (4ttD<R>) 3 V -2 (^) 3 / dk , (k') 2 (-1/D(k' ) 2 ) 2 / d* k , / dO^sine^ 

o o o 

oo 

= ( 4ttD<R> ) 3 V~ 2 • . / dk'fk')” 2 = » . (2.3.24) 

(2tt) D o 

The integral on the third line of Eq. (2, 3.24) is elementary/ and also 
infinite ! We deal with this difficulty below* 

The left-hand sides of Eqs. (2*3*18) and (2.3*24) may be repre- 

A 

sented in terms of diagrams* Since G(k) is the Fourier transform of the 
Green's function for free diffusion (which occurs between the grains)/ we 

A 

represent a factor of G(k) by a line containing an arrowhead. On the other 

A 

hand/ T^(k) is the Fourier transform of the operator describing interaction of 
a solute molecule with a precipitate grain of radius • We represent a 

A 

factor of T^(k) as a circle* Hence, the left-hand side of Eq. (2.3*18) may be 
represented as shown in Figure 2.2a 



i « 


Figure 2*2a. Diagram representing the left-hand side of Eq* (2*3*18). 

This may be interpreted as interaction at grain "i," followed by free diffu- 
sion, followed by interaction at grain "j,” followed by free diffusion, fol- 
lowed by interaction at grain M £* H The left-hand side of Eq* (2*3.24) can be 
represented as shown in Figure 2*2b, 
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Figure 2,2b. Diagram representing the left-hand side of Eq. (2.3.24). 

This may be interpreted as interaction at grain "i ," followed by free diffu- 
sion, followed by interaction at grain "j," followed by free diffusion back to 
"i", followed by interaction at "i. w For three grains, this diagram in Figure 
2.2b is generalized as in Figure 2.2c, 


i 



Figure 2.2c. Ring diagram for three grains. 

Closed diagrams of the sort represented by Figures 2.2b and 2.2c are called 
"ring diagrams." Hence, the terms in the triple sum on the right-hand side of 
Eq. (2.3.3) are of two types: those represented by open diagrams, such as 
Figure 2.2a, and those represented by ring diagrams, such as Figure 2.2b. 

Separating the two types of terms, the triple sum may be written 
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NNN aaaaa 
E E E «T (k)G(k)T.(k)G(k)T(k)>> 
i-1 j-1 £-1 1 3 * 


«T. (k)G(k)T . (k)G(k)T (k)» 
3 


N N A A A A A 

+ E E «T. (k)G(k)T . (k)G(k)T (k)>> . (2.3.25) 

i=1 j=1 1 3 1 

*x 

Making a count of terms of the type represented by Figure 2.2a, 
we note that there are N integral values that can be taken on by i, but once 
it is fixed, only (N-1) values can be assigned to j, and once j is fixed, only 
(N-2) can be assigned to £• Hence there are N(N-1)(N-2) terms in the triple 
sum of this type. Counting terms of the second type, i may be fixed in any of 
N ways, which leaves (N-1) ways of fixing j. Hence, there are N(N-1) terms of 
this type in the triple sum. 

Because of Eq. (2.3.18), each term in the triple sum on the 
right-hand side of Eq. (2.3.25) is identical to every other term in that 

sum. Likewise, because of Eq. (2.3.24), each term in the double sum on the 

right-hand side of Eq. (2.3.25) is identical to every other term and is also 

infinite* Hence, using Eq. (2.3.18), Eq. (2.3.25) may be rewritten 

N N 

Z E «T (k)G(k)T .(k)G(k)T (k)>> 

i-i j-1 1 j * 


N 

E 

i=1 


N 

N 

N 

= E 

E 

E 

i=1 

j-1 

£=1 








= (|K ! v 1 )(^)wrfKR» 3 (S(k)) 


2 


+ 


N(N-1) « T 1 (k)G(k)T 2 (k)G(k)T 1 (k)>> 
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= p 3 ( 4 ttD<R> ) 3 (G(k) ) 2 + N(N-1) « ^ (k)G(k)T 2 (k)G(k)T (k)>> (2.3.26) 


where in the final step, we have taken the thermodynamic limit. In the second 

A 

term on the right in each line of Eq. (2.3.26), we have replaced I\(k) by 
T^ (k) and T_.(k) by T 2 (k), because Eq. (2.3.24) is independent of i and j; 
hence, we may choose any convenient values for i and j so long as they are 
different. 

We now substitute Eq. (2.3.26) into Eq. (2.3.3) and use Eqs . 
(2.3.11), and (2.3.16) to simplify the results. We obtain 


K t (3) + 4ttD p <R>G( k) 4 ttD p <R>G( k) 4 ttD p <R> 

= p 3 (4 ttD<R> ) 3 (G( k) ) 2 + N(N-1)<<T 1 (k)G(k)T 2 (k)G(k)T 1 (k)» . (2.3.27) 


Because the terms involving <R> on either side of Eq. (2.3.27) cancel, we have 

K t ( 3) = N(N-1)<<T 1 (k)G(k)T 2 (k)G(k)T 1 (k)>> (2.3.28) 

which shows that 1^(3) is determined solely by the infinite terms in the 
triple sum. 

2.3.4 Summing Ring Diagrams 

So far, our analysis has carried us only up to the partial rate 
coefficient 1^(3). If we were to push on to find 1^(4), we could do so in 
complete analogy with the above. We would find 


N N N 

K (4) = E E E «T (k)G(k)T.(k)G(k)T (k)G(k)T (k)» . (2.3.29) 

i=1 j=1 1 3 1 1 


The right-hand side of Eq. (2.3.29) is a sum over ring diagrams involving 
three grains. The diagrams are of the type shown in Figure 2.2c. When the 
ring sum in Eq. (2.3.29) is reduced as in the steps leading from Eq. (2.3.25) 
to Eq . (2.3.28), we f i nd 
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K t ( 4) = N(N-1)(N-2)<<T 1 (k)G(k)T 2 (k)G(k)T 3 (k)G(k)T 1 (k)>> 


(2.3.30) 


Now read Eq. (2.3.30) from right to left. Refer respectively to Eqs. (2.2.28) 

A 

and (2,3.21) for the definition of (k) and for an example of the effect of 

the operator P(k|k* ). In reading Eq. (2.3.30) from right to left, the first 

A 

T^k) operates on unity to its right, which has no effect. Hence, this opera- 
tor is just a function of k in no wise different than G(k) . The next operator 

A 

encountered in this direction is T (k) which replaces k by k' in the product, 

A A J A A 

G(k) T^(k) on its right* Once this is accomplished, the product G(k)T^(k) is 
no longer a function of k although T (k) is* The next operator encountered is 

A J A A 

T^(k), which changes k to k" in the product G(k)T^(k) to its right. 

However, k* and k" are independent variables of integration and may be inter- 
changed (i.e., whichever is immaterial). This means that the order in which 

A A A A 

the factors of T^(k)G(k) and T^(k)G(k) occur in Eq. (2.3*30) is unimportant, 
and we may group them in arbitrary fashion. Moreover, the integrals implied 

A A A A 

by the terms T^(k)G(k) and T^(k)G(k) are the same so we may write 

1^(4) = N(N-1)(N-2)<<T 1 (k)G(k)(T 2 (k)G(k)) 2 T 1 (k)>> . (2.3.31) 

The right-hand side of Eq. (2.3.31) is infinite just as was found in Eq. 
(2.3.24) . 

If we were to generalize Eq. (2.3.29) still further to compute 
K^(T ), we would confront rings of the type shown in Figure 2. 2d 



Figure 2. 2d. Ring diagram for four grains, 
and also multiply connected rings of the type shown in Figure 2.2e. 




Figure 2.2e. Multiply connected ring diagram for four grains. 

The diagrams in Figure 2.2e are less singular than those shown in Figure 2. 2d 
and will be ignored in our subsequent analysis. 

Hence, our final result is a sum of terms of the general type 
displayed in Eqs. (2.3.11), (2.3.28), and (2.3.31). That is 


K t = K t (1) + K t (3) + K t (4> + ... 


(2.3.32) 


where Kp(2) is missing because of Eq. (2.3.16). The sum on the right-hand 
side of Eq. (2.3.32) is diagrammed in Figure 2.3. 



Figure 2.3. Diagrammatic representation of the sum on the right-hand side 

of Eq. (2.3.32). 

A A 

If we add and subtract the quantity N<<T^ (k) >>G(k) >> to the right-hand 
side of Eq. (2.3.32) and use Eqs. (2.3.10), (2.3.28) and (2.3.30), we obtain 

= N«T 1 (k)>> + N<<T 1 (kjGfk)^ (k)>> - N<<T 1 (kjGfk)^ (k)>> (2*3.33) 


+ N(N-1)«T 1 (k)G(k)(T 2 (k)G(k))T 1 (k)>> 
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+ N(N-1 ) (N-2)<<T 1 (k)G(k) (T 2 (k)G(k) ) 2 T 1 (k)» 


By Eq. (2.1.5), we imply that K(R) is the effective rate coef- 
ficient for a single grain, the configurations and radii of all other grians 
taking on average values. We may obtain a formula for K(R) at R = R-j simply 
by failing to take the average over P(R 1 )dR 1 implied by the brackets in Eq. 
(2.3.33). Also factoring out the common factor N and grouping other terms, we 
may write Eq. (2.3.33) in the form 

K(R 1 ) = (-^)/d 3 r 1 |T 1 (k) + T^kJGWfl + (N-1 ) («T 2 (k)G(k)» ) 

+ (N-1)(N-2)(«T 2 (k)G(k)>>) 2 + • • • (k) - T. J (k)G(k)T 1 (k)} . (2.3.34) 


Referring to Eqs. (2.3.7) and (2.3.8), <<T 2 (k)G( k) >> may be computed as 
follows: 

«T 2 (k)G(k)>> = 

00 

f dR 2 P(R 2 )47rDR 2 v' 1 f 2 ^) 3 fd 3 k' fd 3 r 2 e' l(k ' -10 ’ r 2 P(k|k’ )G(k) 
o 

oo 

= / dR 2 P(R 2 )4TTDR 2 v' 1 f^:l 3 {d 3 k'(2Tr) 3 6(k , -k)G(k') 

o 


= 4TrDV _1 <R>G(k) . (2.3.35) 

where in going from the first to the second lines, we have as usual inter- 

3 3 

changed the order of integration with respect to d k 1 and d r 2 and recognized 
a delta function. Substitution of Eq. (2.3.35) into Eq. (2.3.34) yields 

K(R 1 ) = + T^kJGfkjri + f 4irD<R>G ( k) ) 
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+ (4nD<R>G(k) ) 2 + ...jT^k) - T^klGtklT^k)} 

= p/d^jT (k) +T (k)G(k)[l + p(4irD<R>G(k) ) + p 2 (4irD<R>G(k) ) 2 


A A A 


+ . . . ]T 1 (k) - T^kjGfkjT^k) } 


(2.3.36) 


where in the last line we have taken the thermodynamic limit* 

In Eq. (2.3.36 ) , the series within the square bracket is a geo- 
metric series in the variable, 4ttDp<R> / and sums to (1 - (4irDp<R>)) ^ • We 
substitute this result into Eq. (2.3.36) and simplify. 


K(R 1 ) = p /a 3 ^ {T t (k)+T 1 (k)G(k) [ 1 - 4irDp<R>G( k) l" 1 T^k)-T 1 (kJGfk) 1 ^ (k) } 


= p/d 3 r 1 {T 1 (k) + T (k)G(k)[- 


1 


1 - 4irDp<R>G(k) 


- 1 1 T , (k) } 


= p/d 3 r {T (k) + T (k)G(k) [ • 4ir ?P <R >?^ ) ]t (k) } 

1 - 4TrDp<R>G(k) 


p/d 3 r T (k) + p/d 3 r T (k)G(k) [ . - i2L D P<S > G < ^> ]t (k) . 

1 - 4irDp<R>G(k) 


(2.3.37) 


The first integral on the right-hand side of the last line in Eq. 
(2.3.37) is the same as Eq. (2.3.8) except for a factor of V -1 and an average 

** 1 °° 

over R^. Removing V fdR^PfR^) fran Eq. (2.3.8), i = 1, and calling the 

o 

resulting integral I^, we have 


I 1 = /d 3 r 1 T 1 (k) = 4 ttDR 1 


(2.3.38) 
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In the second integral on the right-hand side of Eq. (2.3.37), we 

A A 

replace G(k) within the bracket by Eq. (D.35), replace (k) by its definition 
in terms of Eq. (2.2.28) and simplfy. In evaluating the effect of the 
P(k|k' ) operators, we rely upon Eq. (2.3.21). Calling the integral I 2 , we 
find 


i 2 - r«V l( kwk)f no 

1 - 4irDp<R>G(k) 


= /d 3 r 1 (4 1 rDR 1 )(-^-) 3 /d 3 k* 


i(k'-k).r. 
e 1 


ptkik'wwt 4 ’ D p< R ><-vk °> 


1 - 4irDp<R>(-1/k D) 


x (4 ir DR 1 )(-^-) 3 /d 3 k"e l(k "" k) * r 1 P(k|k") 


= ^) 3 /A.a ( k.> 


( 4-irDR^) -4irp<R> ) 
(k' ) 2 + 4irp<R> 

(4ttDR ) 2 (-4irp<R>> 
(k' ) 2 + 4 tt P <R> 


(^) 3 /dV|d 3 r 1 e i,k "- k ’ ,r i 

(2”] 3 /d 3 k" ( 2 it ) 3 (S ( k"-k) 




(4irDR ) 2 (-4irp<R>) 
(k' ) 2 + 4irp<R> 


(2.3.39) 


According to Eq. (D.35) G(k* ) depends only upon the length of k 1 ; hence, the 
integral in Eq. (2.3.39) is spherically symmetrical, and it is especially 
worthwhile to perform the integration in spherical polar coordinates. We find 


' ( 2 I ) 3 / > 2 (^4; 

" O D(k’) 


( 4irDR 1 ) ( -4trp<R> ) 
(k* ) 2 + 4irp<R> 
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2 7T tt 

/ ^ n ^k* 

o o 


( 4ttDR 1 ) (4tt P <R>) 

( 2 TT ) 3 D 


(4n) • / 

o (k' ) + 4irp<R> 


(2.3.40) 


The integral over k' on the last line of Eq. (2.3.40) can be found in tables 

- 

and is equal to ( tt/2 ) (4ttp<R> ) • Hence Eq. (2.3.40) becomes 


(4irDR ) (4irp<R>) (4 tt) 

(2ir ) 3 D 


(tt/2) 


(4ttp<R>) 


1a 

= (4 7 rDR 1 )R 1 (4 T rp<R>) 1 . 


(2.3.41) 


Finally , we substitute our evaluations of I 1 and I 2 [Eqs. 

(2.3.38) and (2.3.41)] into Eq. (2.3.37) to obtain our final result for K( R^ ) , 


which is 


K(R 1 ) = 4ttDR [1 + R 1 (4ttp<R>) /2 ] . 


(2,3.42) 


In the limit of an infinitely dilute system of sinks (p+0), we have according 
to Eq. (2.3.42) 


lim K(R J| ) = 4 ttDR 1 


(2.3.43) 


which is identical to Eq. (1.5.27). On the basis of the arguments and results 
of this section, we have demonstrated that for a slightly concentrated system 


of sinks that corrections to the infinitely dilute limit rate constant. 




v 2 

4ttDR 1 , depend upon p S 

We now return to Eqs. (1.5.26) and (1.5.35), replacing in the 
former the factor 4ttDR by Eq. (2.3.42). The resulting equation is 


v(R ' C A bUlk > ■ 3? ■ (’ + W‘>"P<R»' /2 Uc A bulk - c A (W] 

A 


( 2 %3 %44 ) 


which gives the rate of growth of the grain radius taking into account the 
proximity effects associated with neighboring grains. 

If the bulk concentration is changing with time, i.e., = 

C A (t), Eq. (2.3.44) may be rewritten 

V(R,t) = [1 + R(4irp<R>) 1/2 1 [C (t) - C (R) 1 . (2,3.45) 

x a r 

During a precipitation, R = R(t) will be an implicit function of t. Eq. 
(2.3.45) shows that V(R,t) = V(R(t),t) depends implicitely on t through R(t) 
and explicitly on t through C A ( t) . 
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III. DETERMINATION OF THE GRAIN-SIZE DISTRIBUTION FUNCTION 

AND RELATED PARAMETERS 

3.1 Summary of Governing Equations 

In Section 1.2, we found that the local solubility in the parent phase of 
a component A just outside a grain of radius R was given by 

C (R) - C (oo)exp(2YV;/kTR) • (3.1.1) 

A A A 

The grain size distribution function, n(R(t),t), was introduced in Section 1.3 
and was found to satisfy the continuity equation 


3n(R(t) ,t) 9(n(R(t) ,t)V(R(t) ,t) ) = 

8t 3R 


(3.1.2) 


In Section 1.4, the equation. 


Cl (0) = C (t) + (— ) X* / dRR 3 n(R,t) 
A A v . A J 


(3.1.3) 


3V 1 


was derived expressing the conservation of component A in the case where the 
grains were concentrating A while the parent phase was rejecting it. Finally, 
at the end of Chapter II, we found that the grain radius growth velocity was 
given by (see Eq. (2.3.45)) 


V(R( t) , C ft (t)) = ^ [1 + R(4rrp<R>) 1/2 1 [cyt) - C A (R)] 
A 


(3.1.4) 


where p is the number density of grains at time t. 

The object of the remainder of this Chapter is to solve Eq. (3.1.2) for 
n(R(t),t) subject to the constraint provided by the mass conservation 
principle expressed by Eq. (3.1.3). To start the process of solution, Eq. 
(3.1.1) will be used to evaluate C A (R) in Eq. (3.1.4) and the resulting 
formula will be substituted for V(R(t),t) in Eq. (3.1.2). 


52 


3.2 Introduction of Dimensionless Variables 


Since the parent phase is rejecting the substance A, the concentration 
C^(t) is a monotonically decreasing function of time. When t-n», the 
precipitation is complete, and the parent phase has reached the saturation 
concentration Prior to that, C^( t) > C^foo), which permits us to define 

a diminsionless supersaturation, a(t), by 


cr{ t) 


C A (t> - C A ( ”> 

V“> 


(3.2.1) 


By the time, t+co, at least one grain has achieved macroscopic dimensions? 
mathematically expressed, this means The equilibrium solubility of a 

macroscopic grain is C^(oo) given by Eq. (3.1.1), which must be identical to 
the long time limit of C A (t). Hence, 


lim C (t) = lim C (R) = C («) . (3.2.2) 

A A A 

t*><» R+co 


We conclude that symbols C (<») in Eqs . (3.1.1) and (3.2.1) are identical. 

The number of grains per unit volume of the system with radii between R 
and R + dR is n(R,t)dR. Hence, at any time the grain number density is 


00 

p(t) = / dR n(R,t) . (3.2.3) 

o 

3 

The volume of a grain of radius R is (4 tt/ 3)R . Hence, the fraction of 

00 

the total volume occupied by grains is f dR(47r/3)R n(R,t). As t-n», we may 

define a limiting volume fraction A by° 

00 


<J> = (4tt/3) / dR R 3 n( R,<») = (C (0) - C _ f «>) ) ( V' /X • ) 

' oo J A A A ’ 

O 


(3.2.4) 


where to obtain the far right-hand side of Eq. (3.2.4), we have used the limit 
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of Eq. (3*1.3) as t->»«>. 

We now define a dimensionless length, a, by 


a = R/a 


(3.2.5) 


where 


a = 2 Y V^/kT 


(3.2.6) 


and a dimensionless time, x* by 


t = tD(VVX^)C A (oo)/ci . 


(3.2.7) 


We convert from n(R,t) to a dimensionless distribution, F(a,x), by writing 


n(R,t) = ( 


4ira^/3 




(3.2.8) 


Substituting Eqs . (3.2.1), (3.2.4), and (3.2.7) into Eq. (3.1.3) gives 


cr(x) - a(0) [1 - / da F(a,x)a ] 

o 


(3.2.9) 


where 


V 0) - t V“> 


Turning to Eq. (3.2.3), we substitute Eq. (3.2.8) and obtain 


p(t) = ( ^ — ) p( t) 

4ttcx /3 


where 


(3.2.11) 
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00 


(3.2.12) 


p(x) = / daF( a,x) . 
o 

It is helpful also to define 
00 

/ daF(a,x)a 

<a> = — . (3.2.13) 

00 

/ daF( a, t) 
o 

Next in Eq. ( 3 • 1 • 1 ) r note that 2yV^/kTR is much less than unity for most 
substances at room temperature in the form of precipitate grains exceeding a 
few hundred Angstroms in radius* Hence , we may expand the exponential in Eq. 
(3.1.1) and write more simply 

C (R) = C (oo) (1 + ot/R) . (3.2.14) 

A A 

Substitute Eq. (3.2.14) into (3.1.4) to obtain 

V( R ,t) = M + R(4 tt P <R>) 1/2 ] [(C A (t) - C A (-)) - aC A (oo)/R] . (3.2.15) 

A 

Next use <R> = a <a>, replace p by Eq. (3.2.11), and replace C ft ( t) - C A (°o) by 
Eq. (3.2.1). With these substitutions, Eq. (3.2.15) becomes 

? ' D C A (W) ~ 1/2 

V(a, T ) = -ri — (a(t) - (1/a)] [1 + a(3«J, p<a>) ' ] . (3.2.16) 

a a °° 

2 - 

Next substitute dR = ada, dt = X'ct dt/DV'C («) , Eq. (3.2.8) and Eq. 

A A 

(3.2.16) into Eq. (3.1.2). The result is 

+ |x {F(a,x)(|)[a(T) - + a(3^p<a» 1/2 ]| = 0 (3.2.17) 
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which may also be written 


3F(a,x) + 3 [F(a / T)W(a / T)l _ Q 
3t 3a 


(3*2*18) 


where W(a,x) is a dimensionless velocity defined by 


W(a,T) = (^)fa(x) - (“)][1 + 3< t> 00 P <a> ) 1 • 


(3.2.19) 


Summarizing the accomplishments of this section, we may state that we 
have combined Eqs* (3*1*1) and (3*1.4) and substituted the result into Eq. 
(3.1.2) which has been expressed in dimensionless variables to obtain Eq. 

(3.2.17) . We have also written Eq. (3.1.3) in dimensionless variables to 
obtain Eq. (3.2.9). The four equations (Eq. (3* 1*1)— (3.1.4) ) have been 
condensed to two equations (Eq. (3.2.9) and (3.2.17)). We must now solve Eq. 

(3.2.17) for F(a,x) subject to the side condition expressed by Eq. (3.2.9). 

3.3 Similarity Transformation 

Eq. (3.2.18) is a rather complicated first order partial differential 
equation, whose exact solution does not seem to be known. Hence, we propose 
to separate variables and search for a series solution of the form good for 
large t, 

-y 

F ( a, t) = T Y [F (Z) + x F ( Z) + ••••] , (3.3.1) 

o 1 


where y-j>y>0. The variable Z is defined by the similarity transformation 

a( t) = Z( t )t X (3.3.2) 

where x is some power. If for large t/ where Eq. (3.3.1) is meant to apply, 

55 

we suppose that a is growing as a power of T/ say a ~ t t then Z can depend 
only weakly upon x« In this sense, it may be said to be the proper or 
canonical length for the problem. 
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Substitute Eq. (3.3.1) and (3.3.2) into Eq. (3.2.9). Noting that 
da = r X dZ. The result is 


oo — y oo 

ct(t) = cr( o ) [ 1 - t 4X_Y ( f Z2 3 F q (Z) + T 1 f dZZ 3 F 1 (Z) + (3.3.3) 


' y 1 

For large x, o(x) ~ 0. Since t does to zero as t+<», Eq. (3.3.3) can be 
satisfied for large t only if 


y = 4x 


(3.3.4) 


and 


f dZZ 3 F (Z) = 1 . (3.3.5) 

J o 

o 

Equation (3.3.5) is required in order to cancel the leading term of unity 
within the brackets in Eq. (3.3.3). Using Eqs. (3.3.4) and Eq. (3.3.5)/ Eq. 
(3.3.3) may be written to lowest order as 


a(x) 



(3.3.6) 


where 


a 1 = -a(0) / dzz 3 F 1 ( Z) . (3.3.7) 

o 

Keeping only the leading terra in Eq. (3.3.1) and using Eq. (3.3.4), we 
obtain 

F(a,x) - t~ 4x F (Z) . (3.3.8) 

o 

Substitute Eq. (3.3.8) into Eq. (3.2.12). The result is 


p( t) 


s P, 


-3 x 
T 


(3.3.9) 
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where 


03 

p = / dZ F (Z) # 
K o J o 

o 


Likewise substitute Eq* (3.3.8) into Eq. (3.2.13) to obtain 

. v x 

<a> * a x 
o 


where 


a 

o 


oo 

f dZZ F (Z) 
J o 

o 

oo 

/ dZ F £Z) 
o 


We next introduce the canonical length into Eq. (3.2.18). 
the relevant derivatives, however, it is important to note that 
implicit function of t* Using Eq. (3.3.2), we compute 


9F( a ,x ) 
3t 


= — (- 
3t 1 


-4x 


F o (Z) ) 


= -4x 


- (4x+1 ) 

T 


F (Z) + x 
o 


-4x 


3F (Z) 
o 

3Z 


3Z 

3t 


According to Eq. (3.3.2) 


3Z 

3t 



-(x+1) 
-xt a 


-1 

-XT X 


so that Eq. (3.3.13) becomes 


3F(a, T ) 

3t 


-4x t 


- ( 4 x+ 1 ) 


F (Z) 
o 


xZt 


- ( 4 x+ 1 ) 


3F ( Z) 
o 

3Z 


(3.3.10) 


(3.3.11) 


(3.3.12) 

In computing 
Z is an 


(3.3.13) 


(3.3.14) 


(3.3.15) 
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The next derivative to compute is 


a[F(a,-r)W(a,T)] = 9[F(Z,t)W(Z, T )1 9Z = -4x 9 [F q ( Z) W( Z,t) ] 9Z 
3a 3Z 3a 3 Z 3a 


(3. 


According to Eq. (3*3*2), 


3Z _ -x 
3a T 


(3. 


so that Eq. (3*3*16) becomes 


MFU.T.W.a.TH -5x 

3a T 3Z 


(3. 


Substituting Eqs. (3.3.15) and (3.3.18) into Eq. (3.2.18), we obtain 


,, 3F (Z) _ 3 F (Z)W(Z, T ) 

~ T [4xF (Z) + xZ — — I = t — 


3Z 


(3. 


which suffices to transform Eq. (3.2.18) to canonical variables. 

The explicit form of the velocity which appears in Eq. (3.3.19) may 
obtained by substituting Eqs. (3.3.2), (3.3.9), and (3.3.11) into Eq. 
(3.2.19). The result is 

-1 1 /2 

W(Z, T ) = (Z T ) [a( t) - (1/Z T )] [1 + (3<f> a ) 7 Z] . (3, 

T 00 'o O 

If we set 


e = (3<J> p a ) 
T oo O O 


fa 


(3, 


Eq. (3; 3.20) may be written simply 


W(z,x) = (1/Z T X ) [<x( t) - (1/Zx X )] [1 + eZ] . 


(3, 


3.16) 


3.17) 


3. 18) 


3.19) 


be 


3.20) 


3.21) 


3.22) 
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The asumptotic form of a(x) is given by Eq. (3.3.6). In Appendix E, we 
showed that = x [see Eq. (E.1.16)]. Combining these results and 
substituting into Eq. (3.3.22)/ we obtain 

W(Z,T) * T~ 2x (1/Z) [a., - (1/Z)] [1 + e z] . (3.3.23) 

Next, we substitute Eq. (3.3.23) into the right-hand side of Eq. (3.3.19). 

_ 2 v 

The factor t which leads Eq. (3.3.23) is combined with the factor 
c — 7x 

of t” :>X on the right-hand side of Eq. (3.3.19) to produce t . Because of 

— ( 4 x+ 1 ) 

Eq. (E.2.2.), which specifies x = 1/3, the factor of x on the left-hand 

-5 X 

side of Eq. (3.3.19) cancels the factor of x on the right. The end result 
of these cancellations is shown in Eq. (3.3.24). 


dF (Z) h -1 

-4xF (Z) - xZ — 2_ — + — [F (Z) Z (a. - (1/Z)) (1 + ez)] 
o dZ dz o 1 


0 . (3.3.24) 


Equation (3.3.24) may be cast in a more compact form by adding -xZF Q (Z) to the 
terms within the brackets and subtracting -(d/dZ) ( xZF q ( Z) from the terms 
outside. The result is 


-F (Z) + — [w(Z)F ( Z) ] = 0 
o dZ o 


(3.3.25) 


where we may identify a new dimensionless velocity by 

»<« - (?)(*, -i)« -fl • 


(3.3.26) 


We now summarize the progress to this point. In Section 3.2, we were 
left with Eq. (3.2.18) to be solved for F( a , t) subject to the condition 
on cj(t) expressed by Eq. (3.2.9). By expanding the grain size distribution 
function as in Eq. (3.3.1) and introducing the variable change specified by 
Eq. (3.3.2), we have converted Eq. (3.2.9) into Eq. (3.3.6) with o ^ defined by 
Eq. (3.3.7). The variable change specified by Eq. (3.3.2) plus the simple 
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asumptotic forms given by Eqs. (3.3.6), (3.3.8), (3.3.9), and (3.3.11) have 
served to convert the partial differential equation, Eq. (3.2.18), to the 
ordinary differential equation, Eq. (3.3.25), with w(Z) defined by Eq. 

(3.3.26). The constraint expressed by Eq. (3.2.9) has been completely 
incorporated into Eq. (3.3.26) and no further refernce to it is needed to 
obtain the general solution of Eq. (3.3.25) for F Q (Z). 

In the process of changing variables from a(t) to z(x)r we have converted 
Eq. (3.2.12) to the form of Eqs. (3.3.9) and (3.3.10) and converted Eq. 
(3.2.13) to the form of Eqs. (3.3.11) and (3.3.12). Once we have solved Eq. 
(3.3.25) for F q (Z), we will use Eq. (3.3.10) to compute the normalization 
factor for F q (z) and use Eq. (3.3.12) to compute the mean radius. The 
technique used to solve Eq. (3.3.25) is motivated by some physical arguments 
which we explore in the next section. 

3.4 Critical Radius and Supersaturation 

As was explained in Section 1.1, in a distribution of grain sizes 
maturing through Ostwald ripening, there exists a critical radius specified by 
the value of the super saturation. Grains smaller than this critical radius 
will dissolve while larger grains will grow. As the precipitation evolves, 
however, the supersaturation diminishes and the critical radius increases. In 
the case of a grain whose initial radius exceeds the critical radius, there 
ensues a race between it and the critical size. Even though such a grain is 
initially growing, it is possible for the critical radius to grow faster and 
at some time to catch up. When this happens, a grain of this size stops 
growing, turns around and begins dissolving. In fact, when x=«>, all grains 
except the largest have disappeared. The sole remaining grain has absorbed 
all of the available solute and is in stable equilibrium at the saturation 
concentration. 

For large yet finite times (i < <») there is an asymptotic distribution of 
grain sizes given by the function F q (Z). Within this distribution, a grain 
which is neither growing nor dissolving has attained the critical radius. 
Because it is in equilibrium with the instantaneous superstauration, its 
growth velocity, as defined in canonical variables by Eq. (3.3.26) must be 
zero. The canonical critical radius Z Q thus satisfies 

w(Z Q ) = 0 « (3.4.1) 
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Combining Eqs. (3,3.26) and (E.2.2) and substituting the result into Eq. 
(3.4.1) , we obtain 

(1/z )(a- - (1/z ))(1 + eZ ) - (1/3) Z = 0 . (3.4.2) 

o 1 o o o 

Equation (3.4.2) may be rearranged to read 


a 


1 



o 


3( 1+eZ ) 
o 


+ 


Z 


(3.4.3) 


so as to serve as an expression for a ^ • 

Now consider a grain with radius , Z > Z Q . This grain cannot grow forever 
or else the principal of mass conservation expressed by Eq. (3.3.5) would be 
violated, and no asymptotic steady state could exist. Hence, asymptotically 
no velocity larger than the critical velocity, w(Z Q ) = 0, can exist. 

Therefore w(Z) must have a maximum at Z Q . The condition specifying this is 

= 0 at Z = Z . (3.4.4) 

dZ o 

Using Eq. (E.2.2) and substituting Eq. (3.3.26) into Eq. (3.4.4), we 
obtain 




(3.4.5) 


Equation (3.4.5) may be solved so as to provide another expression for o ^ in 
terms of Z Q , namely. 



o 


(3.4.6) 
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Eliminating a ^ between Eqs. (3.4.3) and (3.4.6) yields a polynomial 
equation for Z Q , 

eZ 4 + 2Z 3 - 3e 2 Z 2 - 6 e Z -3=0. (3.4.7) 

o o o o 

For e = 0, Eq. (3.4.7) becomes a simple cubic equation, the physically 
significant root of which is 

Z Q = (3/2) 1/3 . (3.4.8) 

Equation (3.4.8) serves as a zeroth order estimate of the root of the quartic 

3 

which is Eq. (3.4.7). Isolating the term in Z q on the left-hand side, Eq. 
(3.4.7) may be rewritten 


Z 

o 



[' 


( 22 o - -T-) 


e + z. 



(3.4.9) 


Now substituting Eq. (3.4.8) for Z Q on the right-hand side of Eq. (3.4.9), one 
obtains 



(3.4.10) 


Expanding the right-hand side of Eq. (3.4.10) in a binomial series, including 
terms through order z, yields 

, V3 , , V3 

Z o*(f) t 1 + *!• 


Replacing e by Eq. (3.3.21), Eq. (3.4.11) becomes 


Z 

o 



1 - 1/3 1/2 

t 1 + l(f) (3 ^ p o V ] • 


(3.4.12) 
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To calculate substitute Eq. (3.4.11) for Z Q in Eq. (3.4.6)* The 

result is 


9 1/3 i o V3 -1 1 ~ 2/3 , - 1/3 2 

”, - 2 (|) [’ - (|l(|) el - [' e] + e • <3.4.13) 


The quantities within brackets are next expanded in bimonial series up to 
order e. When coefficients of like terms are collected, one obtains 


, V3 . , V3 

0, - (|l [' - Cjlffl el . 


(3.4.14) 


Replacement of e by Eq. (3.3.21) gives finally 


9 1/3 . 4 V3 1/2 

”,*(!) <3*.p 0 . 0 > ,/a i • 


(3.4.15) 


3.5 Integration of Eg. (3.3.25) to obtain F q (Z) 

In Eq. (3.3.25) the indicated differentiation of the quantity within the 
brackets may be carried out and the result divided through by w(Z)F q (Z). One 
thus obtains 


, JF o< 21 

F (Z) dZ 
o 


_1 dw(Z) _J 

w(Z) dZ + w(Z) 


(3.5.1) 


Integration of Eq. (3.5.1) yields 


F (Z) 
o 


C 

o 


w(Z) 


exp [ / 


dZ 

w(Z) 1 


(3.5.2) 
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where C Q is a constant of integration* 

To compute the integral indicated in Eq. (3*5*2), we begin by combining 
Eqs* (3.3.26) and (E.2.2) to obtain 

w(Z) = (£)[ CTi - (£)] (1 + eZ) . (3.5.3) 

For on the right-hand side of Eq* (3.5.3), substitute Eq. (3.4.6). The 

result is 

z 2 

w(Z) = (I)[| | + e - ~1 (1 + eZ) - | . (3.5.4) 

o 

Multiply the right-hand side of Eq. (3.5.4) by (-3Z 0 Z 2 )/(-3Z Q Z 2 ) to obtain 


w(Z) = f -" 1 -■ -) fzf(Z 3 - 6) + 3(Z /Z) - 3Z e I f 1 + eZ] + Z Z 3 } . (3.5.5) 

v ^ /LL o o o JL J o J 

3Z Z 
O 

If on the right-hand side of Eq. (3.5.5) terms of order e and higher are 
dropped, the result is 

w(Z) a [Z z 3 + (Z 3 - 6) + 3Z Q ] . (3.5.6) 

3Z Z ° 

o 

If terras of order e and higher are ignored in Eq. (3.4.10), Eq. (3.4.8) 

3 

applies. Since according to Eq. (3.4.8) 3 - 2Z = 0, we may add 

3 2 ° 

(3 - 2Z^)Z to the quantity within the bracket in Eq. (3.5.6) to obtain 

w(Z) ^ f ) fz Z 3 + (3 - 2Z 2 ) Z 2 + (Z 3 - 6) Z + 3Z ]. (3.5.7) 

v 2 ' L o o o o 

3Z Z 
o 

The quantity within the bracket in Eq. (3.5.7) factors so that 
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(3.5.8) 


w(Z) a -(1/3Z 2 ) (Z - Z ) 2 (Z + 3/Z 2 ) . 

o o 


Evaluation of Eq. (3*5.2) calls for 1/w(Z). Also involved is an integral of 
this quantity. To compute the integral, we expand 1/w(Z) in partial 
faractions and arrive at 


3Z 


w(Z) 


(Z 3 +3) 

o 


(z-z ) 

o 


3z 3 (z 3 46) 
o o 

3 2 

(Z J +3r 

o 


27 


(Z-Z o ) 


3 2 

(Z +3) 
o 


( Z+3/Z ) 
o 


(3.5.9) 


The details associated with the derivation of Eq. (3.5.9) may be found in 
Appendix F. 

After substitution of Eq. (3.5.9) into the integral on the right-hand 
side of Eq. (3.5.2) and integrating, one obtains 


f dz 
J w(Z) 


4 3 3 

3Z „ 3Z (Z* 5 +6) 

3° * j ° ~ 3 P ' ' 2 * n(Z 0 _Z) 

(Z 3 +3) (Z o Z) (Z 3 +3) 2 ° 

o o 


~ — - Jln(Z+3/Z 2 ) 


(3.5. 10) 


(Z +3) 
o 


In Eq. (3.5.10), we have used the following: 

f dZ(Z - Z )“ 2 = (Z - Z) -1 
J o o 

f dZ( Z - Z ) -1 = fcn(z - Z) 

J O o 

2 -1 2 

/ dZ(Z + 3/Z ) = Zn(Z + 3/Z ) . 

J o o 

Finally, we substitute Eqs. (3.5.9) and (3.5.10) into Eq. 
obtain 


(3.5.11) 

(3.5.12) 

(3.5.13) 


(3.5.2) to 


F (Z) 


C Z 


2 


o 

(Z -Z) a (Z+3/Z 2 )® 
o o 


exp(- 6/(Z -Z)) 
o 


(3.5.14) 
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where a, 8r and 6 are defined by 


3Z 3 Z 3 +6) 

a = ° ° -r + 2 (3.5.15) 

(Z +3) 
o 


e = + 1 (3.5.16) 

(Z 3 +3) 

O 


6 = ? - . (3.5.17) 

(Z 3 +3) 

o 


For later comparison, it is useful to have Eqs. (3.5.15)-(3.5.17) in the 
limit e*0 (infinitely dilute system of grains). In this limit, Z Q is given by 
Eq. (3*4.8) and according to Eqs. (3. 5. 1 5 ) — ( 3 • 5. 17) , 


a = 


U 

3 


(3.5.18) 



(3.5. 19) 


6 = 



(3.5.20) 


U s l ng Eq . (3.5.14), 


lim F (Z) 
o 

e+0 


C 1 Z 
o 


[ ( 3/2 ) 1 /3 -Z] 1 1 /3 [Z+2 ( 3/2 ) 1 /3 ] 7/3 


exp(- 


(3/2) 


1/3 


1 /3 

[(3/2) ,/J -Z] 


■) (3.5.21) 


where C' is the e = 0 limit of C . 

o ° 

One should note that Z = Z Q in Eq. (3.5.14) and Z = (3/2) in Eq. 
(3.5.21) cannot be exceeded by Z, otherwise F Q (Z) would take on complex values 
which are unphysical. Hence Z = Z Q serves as a cut-off in the former while Z 
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= (3/2) 1 ^ 3 is a cut-off for the latter. For values of Z larger than these 
cut-off points, the respective distributions are zero. 

3,6 Evaluation of C' , p , and a for e = 0 
o lo 0 . 

3.6.1 Evaluation of C* 

a 

Equation (3.3.5) may be used to determine C* • We note, however, 

o 

that because of the cut-off at Z = Z' q - (3/2 ) 1/ ' 3 , integration beyond this 
point is unnecessary and Eq. (3.3.5) can be rewritten 


Z' 

° 3 

1 = / dZ Z F q (Z) . (3.6.1) 

o 


If we substitute Eq. (3.5.21) into Eq. (3.6.1) and let Z = (2/3) 2 / 3 x , Eq. 
(3.6.1) becomes 


C' 3 j 


27/8 


dx 


( (3/2)-X 1/3 ) (3+x 1/3 ) 7/3 


exp(- 


1 _ 

1- ( 2/3) x 


1/3 


)• 


(3.6.2) 


The integral on the right-hand side of Eq. (3.6.2) has been evaluated 
numerically by Lifshitz and Slyozov [3] and found to have the value 0.0488. 
Equation (3.6.2) implies 

cr = 3/(0.0488) . (3.6.3) 

Equation (3.5.21) may be transformed to another form by 
multiplying the right-hand side by e and dividing ^ by e. Eq. (3.6.3) 
becomes 


C* 

o 


/e 


3 

( 0 • 0488) e 


22.6 


(3.6.4) 


and Eq. (3.5.21) becomes 
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F (Z) = 
o 


22.6 Z 


exp(- 


[ (3/2) V3 -Z] 11/3 [Z+2(3/2) 1/3 ] 7/3 ^ [(3/2) 1/3 -Z] 


). (3.6.5) 


3*6.2 Evaluation of p 


n 


The quantity can be computed using Eq* (3*3*10) 
account the cut-off at Z = Z^ = (3/2) Eq* (3*3.10) becomes 


Taking into 


p o = 


Z' 

O 

f dZ F (Z) 
J o 

o 


(3.6.6) 


Replacing (3/2)^/^ in Eq. (3.5.21) by the sumbol, Z 1 and substituting the 
result into Eq. (3.6.6), one obtains 


p o = 


z' 

O 

J 

O 


dZ 


C' Z 
o 


Z' 

o 


(Z* 

o 


z> , 1 / 3 <z«zm 7/3 

o 


expf 2 1 


(3.6.7) 


(Z' - Z) 
o 


Now in Eq. (3.6.7), we express binomial terms involving fractional exponents 
as exponential functions of logarithms. Last, we multiply and divide the 
integrand by (-3). The result is 


p o = 


c z ; 

-f~J ^ 


(— 3Z ) 


Z* 


exp[-( — 2-) - (2-)£n(Z* -Z) - f^)iln(Z+2Z )] 


---r- 1 v j l-i 

(Z'-Z) (Z+2Z‘) Z'-Z 

o o o 


z' 


z’ 


“ “ ™f“ f ® % ( ex p[“( - (f)£n(z;-Z) - rfnn(Z + 2 z;)l} 

o Z'-Z 

o 


C 1 Z * z 1 

-§- {exp[-r— 2 -) - r|)£n(z;-Z) - r|l£n(Z+2zyn ° 

Z'-Z o 

o 
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(3.6.8) 


(C'/e) 

o 


9(2) 


1/3 


p = 1.994 
'o 

where in the last line, we have used Eq. (3.6.4). 

3.6.3 Evaluation of a 

The quantity, a Q , can be evaluated directly in terms of its 
definition expressed by Eq. (3.3.12) or indirectly by exploiting the results 
of Section E. 1 of Appendix E. We choose the latter procedure because it is 
simpler. 

On the basis of Eq. (E.1.14), we have letting t-m». 


2 _ 

lim f dZ Z F (Z) V(Z) 
J o 

T-M» o 


“Yi 


lim x 


0 


(3.6.9) 


since y 1 = x = 1/3. Into the left-hand side of Eq. (3.6.9), substitute Eq. 
(E.1.15) to obtain 


/ dz (c^z - 1 + ec 1 z 2 - e z) F (z) = 0 . (3.6.10) 

O 


Since by definition, the average of Z k is 


<z k > 


f 

o 


dZ Z k F (Z) 
o 


(3.6.11) 
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Eq. (3.6.10) becomes 


a 1 <z> 


- 1 


+ ea t < z > 


- e<z> = 0 


(3.6.12) 


Because of Eq. (3.3. 12) , <Z> = a Q . Hence when e 388 0 , Eq. (3.6.12) gives 


a o = Va 1 = (4/9) 1/3 

where we have used Eq. (3.4.14) with e = 0. 


(3.6.13) 


3.7 Evaluation of C_, p , and a for e*0 
— - ■ ■ — g — 

When e*0 r C Q must be evaluated by numerical iteration. The independent 

variable is <f> , which is the volume fraction under consideration. One then 
00 

estimates p Q and a Q by choosing the ^ = 0 values given by Eqs. (3.6.8) and 
(3.6.13), respectively. Having specified <)> , p , and a_, Eq. (3.4.12) gives 
Z Q , which in turn determines a, 6, 6 by virtue of Eqs. (3.5. 15) — (3.5. 17) . At 
this point C Q can be determined from Eq. (3.3.5). This completes the first 
iteration. 

With (j)^ fixed as before, the second iteration is begun by using the 
function F Q (Z) obtained from the first iteration and computing p and a Q using 
Eqs. (3.3.10) and (3.3.12). Next, Z Q is recomputed using Eq. (3.4.12). With 
Z Q specified, C Q is determined as before, which completes the second 
iteration. 

The procedure described above converges nicely after a few iterations. 
Part of the results obtained may also be expressed analytically. Let p Q and 
a Q be given by the zero volume fraction results, Eqs. (3.6.8) and (3.6.13). 
Substituting these into Eq. (3.3.21) gives 


e = 2. 1 34 (,J>J 1/2 . 


(3.7.1) 
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When Eq. (3.7.1) is substituted into Eq. (3.4.11), we obtain 


Z = ( 3/2 ) 1 /3 [1 + 1 . 222 ( <f> ) 1/2 ] . 

O oo 


(3.7.2) 


Likewise , substitution of Eq. (3.7.1) into Eq. (3.4.14) gives 

* (9/4 ) 1/3 [ 1 - 0 . 8 1 5 ( <(> ) 1/2 ] . (3.7.3) 

1 T oo 

To obtain a Q , we turn to Eq. (3.6.12). We define a 2 by 

oo 

f dz z 2 f (Z) 

° 

a = <Z > = — = 0.609 . (3.7.4) 

Z oo 

f dZ F (Z) 
o 
o 

The value, 0.609, on the right-hand side of Eq. (3.7.4) was found by replacing 
the integral in the denominator by Eq. (3.6.8) and computing the integral in 
the numerator numerically using Eq. (3.6.5). With a 2 thus found, Eq. (3.6.12) 
may be rewritten 


a 

o 


1 " ea i a 2 



(3.7.5) 


Into the right-hand side of Eq. (3.7.5), we substitute Eq. (3.4.14). Up to 
terms of order e, the result may be written 


1-a 2 (3/2) 2/3 e 

° (3/2) 2/3 -(3/2)e 


(3.7.6) 


Expanding the denominator of Eq. (3.7.6) in a binomial series and collecting 
terms up to order e gives 
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1/3 


a - i ^ 
o ^9 


(f) [' ♦ ( 2 — 2?3 - « 2 (3/2> 2/3 ) e] 


(3.7.7) 


2(3/2) 


We now substitute Eqs. (3.7.1) and (3.7.4) into Eq. (3.7.7). The result is 


1/3 


a o “ (?) / [1 + 0.739 (^) 1/2 1 


(3.7.8) 


This completes the determination of the relevant parameters in the case e^O. 

3.8 Summary and Discussion 

3.8.1 Results in Terms of Laboratory Parameters 

The use of dimensionless variables for mathematical analysis has 
led us to a number of important results. In this section, we summarize a set 
of parameters (mostly dimensional) which are appropriate for expressing these 
results in a form convenient for laboratory use. 

The mean molecular volume in the A-rich phase (the phase 
consisting of grains which are concentrating A) is 


V* 


X A (MW, » * x b (mw, b 


p'N 


o 


(3.8.1) 


where (MW) A and (MW) g are the molecular weights of A and B, respectively, and 
X ^ and are the mole fractions of these components in the phase. The mass 

density of the phase is p 1 , while is Avagadro* s number. The volume 
occupied by a molecule of A in this phase is 

V A - x^ V’ . (3.8.2) 

Using Eq. (3.8.2), we may calculate the natural length scale 

ot = 2y V^/kT (3.8.3) 

which was introduced in Eq. (3.2.6); a is the radius beyond which the 
curvature of a grain no longer greatly enhances its solubility. 
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The equilibrium concentration (molecular number density) of A in 
the B~rich phase (the continuous phase which is rejecting A) is 


P N 

A K o 


C ( 00 ) — 

A X (MW) + X (MW) 


(3.8,4) 


where p is the density of the phase and and X 0 are the respective mole 
fractions • 

With the quantities defined so far, we may also introduce a 
na tur al time sc ale , 



DV’C (oo) 
A 


2 

Having the units of seconds,- 0 is proportional to the time, a /D, required by 
a molecule to diffuse a distance, a* 

By combining Eqs . (3.2.4) and (3.2.10), we obtain an expression 
for the initial supersaturation, a(0), in terms of the volume 
fraction, 6 , occupied by the A-rich phase at equilibrium, 

CO 


X'* m 

cr (0 ) = . (3.8.6) 

C (oo)V' 

A 


According to Eqs. (3.2.7), (3.3.6), (3.8.5), and (E.2.2), the supersaturation 
decays with time as 


a(t) = a^(0/t)^ 3 


(3.8.7) 


where c^ may be determined from the graph in Figure 3.1 or calculated from Eq 
(3.7.3). 



Figure 3.1, Grain size distribution normalization constant/ p [Eq. 

(3.3.10)]/ supersaturation constant/ a ^ [Eq. (3.7.3)]/ and 
mean radius constant/ a Q [Eq. (3.7.8)]/ shown as functions 
of the equilibrium volume fraction occupied by grains/ <J> 

[Eq. (3.2.4)]. 

According to Eqs. (3.2.4)/ (3.2.7)/ (3.2.10)/ (3.2.11)/ (3.3.9)/ 
and (E.2.2)/ the total number of grains of the A-rich phase per unit volume 
decreases with time according to 


o q(0) 1 

4iraD t 


(3.8.8) 




where p may be determined frcm Figure 3.1. 
o 

As we know, within the population/ p(t) r there is a distribution 
of grain sizes. This distribution/ n(R/t) / depends upon both radius and 
time. Combining Eqs. (3.2.4)/ (3.2.5)/ (3.2.8)/ (3.2.10) f (3.3.8)/ and 
(E.2.2)'/ we may write 


n(R/t) 


q(0) 

(D0/ a 2 ) 


F (Z) 
o 

(4ira^/3 ) 



(3.8.9) 
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where F q (Z) is given by Eq. (3*5.14) and Z is defined by Eqs • (3.2.5) and 
(3.3.2) , namely, 

z = (R/a) ( 0/t) 1/3 . (3.8.10) 

The function F Q (Z) is shown in graphical form in Figure 3.2. Note that with 
increasing volume fraction, F q (Z) decreases in height and broadens. The 

broadening is caused because competition increases the rates of both growth 
and dissolution of grains [see Eq. (2.3.45)]. This means that the large 
grains grow even faster, shifting both the mean and maximum radii to larger 
values. With broadening, however, the height must decrease because F Q ( Z) is 
normalized according to Eq. (3.3.5). Associated with the distribution are a 
number of radii, which we discuss below: 



Figure 3.2. Grain size distribution function, F q (Z). 

There is no sharp minimum radius in the distribution. To 
represent an effective minimum radius, however, we arbitrarily select from 
Figure 3.2, Z = 0.2, the point where F q (z) diminishes to a few percent of its 
maximum value. Setting Z = 0.2 in Eq. (3.8.10), we obtain 
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R . = 0.2 a (t/0) 

min 


(3.8.11) 


1/3 

• 

Inspection of Eq. (3.2.15) reveals that there is a critical value 
of the radius, call it at which V(R,t) = 0. Using Eq. (3.2.1), this 

radius may be written 

R crit = a/a(t) * (3.8.12) 

The value of R cr ^ t follows the inverse of the supersaturation. Grains whose 
radii at time t are greater than R cr jLt are growing while grains whose radii 
are smaller than R cr ^ t are dissolving. 

The mean radius in the distribution is found by combining Eqs. 
(3.2.5), (3.2.7), (E.2.2), (3.3.11), and (3.8.5). The result is 

<R> = a Q a( t/0) 1/3 . (3.8.13) 

The value of a Q depends upon the volume fraction and may be read from Figure 
3.1 or calculated from Eq. (3.7.8). Comparison of Eqs. (3.7.3) and (3.7.8) 
for <p = 0, shows that a (<f> = 0) = = 0) = (4/9)^^. Hence, because 

of Eqs. (3.8.7), (3.8.12), and (3.8.13), we find that <R> = R C rit wlien 
<p M = 0. For ^ > 0, however, Eqs. (3.7.3) and (3.7.8) indicate that R crit > 
<R> • 

The maximum radius in F Q (Z) is Z = Z Q (see Eq. (3.5.14) and the 

discussion following Eq. (3.5.21)]. Setting Z = Z Q in Eq. (3.8.10), we obtain 

R = z a(i;/ 0) 1/3 (3.8.14) 

max o 

where Z Q may be read from Figure 3.3 or calculated using Eq. (3.7.2). In 
general R min < <R> < R crit < R^. 
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Figure 3.3* Maximum grain size constant/ Z Q [Eq. (3.7.2)]/ as a function 
of volume fraction occupied by grains, 4> [Eq. (3.2.4)]. 


Finally, it is worth noting that all quantities which depend upon 

volume fraction (e.g., p(t), q(t), <R>, and F Q (Z)) do so as a 

1 /2 

function of (fl^) (see Chapter II above). This has its origin in the growth 

velocity expressed by Eq. (2.3.45), which encompasses the competition effect. 

1/3 

The square root dependence is in contrast with ( tjj^) which might be expected 
if the competition effect depended upon the mean separation between the 
grains, which is proportional to the cube root of the volume. In the 
treatment supplied above, we have assumed that the surface of each grain is in 
equilibrium with its local environment and that its growth is diffusion- 
controlled. If, however, the relaxation of the grain radius to accommodate 
changes in its local environment is not instantaneous, then the growth of the 
grain is not diffusion-controlled. The volume- fraction effects predicted 
above diminish as the degree of diffusion-control decreases, becoming zero in 
the limit that diffusion is infinitely fast. In this limit, the mean radius 
increases as t 1//2 , the number of grains per unit volume decreases as t 2//3 , and 
the grain size distribution function F Q ( Z) is substantially broadened [5] # 
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3.8.2 application to the System Succinonitrile/Water 

The growth of grains described above involves an evaporation- 
condensation mechanism with molecular transport between grains by diffusion. 

In almost all systems, the continuous phase and the granular phase have 
different densities, however. This means that the evaporation-condensation 
mechanism (Ostwald ripening) may be obscured by the gravitational convection 
and collisional coagulation of the growing grains. Expermentation in the 
acceleration-free environment of an Earth orbiting laboratory removes the 
effect of gravitational convection and permits the free observation of the 
Ostwald ripening of the precipitate grains* 

Phase equilibrium in the two component system, succi nonitrile 
(NCCI^C^CN) and water is summarized in the phase diagram shown in Figure 
3.4. The diagram is marked by a monotectic point (two liquids and solid 
C 4 H 4 N 2 ) in equilibrium, which is characteristic of many two-component metal 
alloys. The system succinonitrile/water has some distinct advantages for 
studying the process of Ostwald ripening. These are: (1) The system is 

transparent to visible light at all compositions and temperatures, so that the 
growth of second phase droplets can be observed photographically. (2) The 
consolute and monotectic points occur near room temperature, so that little 
heat is required to keep the system in a liquid state? in space flight where 
electric power is at a premium, this is an important advantage. (3) The 
thermophysical properties of the system are rather well studied. (4) Near 41 
°C, the droplet and continuous phases have the same density, so that at this 
temperature the Ostwald ripening process can be observed even in Earth* s 
gravity. 

In Table 3.1, we have summarized the thermophysical properties of 
the succinonitrile/water system at 25±1 °C. The primed symbols refer to the 
succinonitrile rich phase (droplets). The coefficient for interdiffusion on 
the water-rich side of the phase diagram (continuous phase) is D. The 
interfacial tension is y. The subscript A refers to succinonitrile while the 
subscript B refers to water. 
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temperature (°C) 


P = 760 mm Hg 



Figure 3.4. Phase diagram for the system succinonitrile/water . 


TABLE 3.1 THERMOPHYSICAL DATA FOR THE SYSTEM C 4 H 4 N 2 /H 2 0 AT 25±1 °C 
A = C 4 H 4 N 2' B = H 2°* Primes on symbols refer to the succinontrile-rich phase. 
Succinonitrile Rich Phase Water Rich Phase 


X* = 0.646 
A 

X* = 0.354 
B 

p 1 — 1.0121 gm/ cm^ 


Other data: y = 2.484 erg/cm^; D = 1.27 x 10*^ 
(MW) Q = 18 gm/mole. 


X A = 0.032 
X B = 0.968 
p — 1.0052 gm/ cm^ 

9 

cm /sec ( MW) A = 80 gm/mole; 
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If the continuous phase is to be water-rich, the initial 
composition of the system must lie to the right in the phase diagram. By the 
lever rule, the initial mole fraction of A, X A °, is related to the equilibrium 
values, X A and X^, by 


n 

n* 


x-- X 


o 

A 



(3.8. 15) 


where n' and n are the number of moles of A-rich and B-rich phases, 
respectively, in the equilibrium mixture. The volumes occupied by the two 
phases at equilibrium are, respectively. 


V' = n'V' 


(3.8. 16) 


and 


V = nV 


(3.8.17) 


where 


V' = (X^(MW) a + x^(mw) b )/ p ’n o 


(3.8.18) 


and 


V = (X a (MW) a + X B (MW) B )/pN Q (3.8.19) 

are the respective average molecular volumes in the two phases, the fraction 
of the total volume occupied by the A-rich phase at equilibrium is 


= v ' = 1 
'« V'+V 1 + (V/V') 


(3.8.20) 
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where by Eqs. (3.8.16) through (3.8.19) 


= £l 


V p 


X A (MW) i +X B (MW, B 

x a <mw) a +x b (mw) b 




o 

A 


X 



(3.8.21) 


The thermophysical data in Table 3.1 are for equilibrium at 25±1 
°C. We specify the desired equilibrium volume fraction of the A-rich phase to 
be <b = 0.1, since this is one of the larger values considered in Figures 3.1- 
3.2 and can be expected to demonstrate the competition effect. Combining the 
data in Table 3.1 with ^ = 0.1 according to Eqs. (3.8.20) and (3.8.21), we 
find that the system should be prepared initially so that its composition is 
X° « 0.055, that is it should be 5.5 mole percent C 4 H 4 N 2 an< * ^4.5 
percent water. This and other results are summarized in Table 3.2. 

From Eqs. (3.8.1), (3.8.4), and (3.8.6) we find that the initial 
super saturation, a(0), is 0.70. This means that the initial succinonitrile 
concentration in the continuous phase is 70% greater than it will be when 
equilibrium is achieved. This immediately dete mines the initial critical 
radius for growth or dissolutin through Eq. (3.8.12) with t set equal to 
zero. We find R cr i t (°) = 1*64 x 10“ 4 ym or 1.64 A. 

The other quantities in Table 3.2 cannot be calculated at t = 0, 
because the formulae, Eqs. (3.8.7), (3.8.8), (3.8.9), (3.8.11), (3.8.13), and 
(3.8.14), are all asymptotic and accurate only when t >> 0. We have chosen to 
evaluate these formulae at a time where the mean radius in the distribution 
achieves 1 ym, which is a size convenient for holographic photography. This 
occurs at t = 89 sec, which may be seen from Table 3.2 vastly to exceed 0 ■= 
30.4 ps. At this time, the minimum radius in the distribution is 0.21 ym, 
while the maximum radius is 1.7 ym. The critical radius, 1.1 ym, slightly 
exceeds the mean radius. Having been initially 70%, the super saturation is 
just 6.8 x 10~ 3 %. There are 2.2 x 10 10 droplets of the C^H^^-rich phase per 
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TABLE 3.2 CALCULATED RESULTS FOR THE SYSTEM C 4 H 4 N 2 /H 2 0 AT 25±1 °C 

AND u =0.1 

oo 

The results in the table are based on the following quantities; the equation 
or figure used to specify each quantity is shown in parentheses next to the 
numerical value: V 1 = 9.52 x 10 23 cm 3 /molecule [Eq. (3.8.1)]. V^* = 6.15 x 

10" 23 cm 3 /molecule [Eq. (3.8.2)]. C^too) = 9.66 x 10 20 molecules/cm 3 [Eq. 

(3.8.4) ]. a = 7.43 x 10 -9 cm [Eq. (3.3.3)]. 0 = 3.04 x 10 -11 sec [Eq. 

(3.8.5) ]. Pq” (Figure 3.1). a 1 = 0.97 (Figure 3.1). a Q = 0.94 (Figure 

3.1). Z Q - 1.6 (Figure 3.3). X° = 0.055. 


t( sec) 

p (cm -3 ) 

a 

*min ( **) 

R C r it <um) 

<R> ( ym) 


0 

— 

7.0 x 10 -1 


1.1 x 10** 4 

— 

— _ 

89 

2.2 x 10 10 

6.8 x 10 -5 

0.21 

1.1 

1.0 

1.7 

3.45 x 10 5 

5.6 x 10 6 

4.3 x 10 -6 

3.3 

17 

16 

27 


After 4 days (345,600 sec) , which might be the length of an 
experiment carried into orbit by the space shuttle, the mean radius has 

increased to 16 \im. The minimum radius is 3.3 yim, while the maximum radius is 

—6 —4 

27 ym. The remaining super saturation is just 4.3x10 =4.3 x 10 %, which 

is extremely close to equilibrium. The number of droplets has decreased to 
5.6 x 10 6 per cm 3 . 

—9 

Because the characteristic length, a = 7.43 x 10 cm, is quite 
small, the approximation, R > a, invoked to obtain Eq. (3.2.14) is established 
early in this system. In fact, because the interfacial energy, y = 2.484 
erg/cm , is small, the approximation becomes good when the droplet has just 
barely achieved molecular size. It is doubtful that interfacial energy even 
can be defined at this small size, although this in no way diminishes the 
power of the approximation. 

Likewise, because a is so small, the characteristic time 
scale, 0 = 30.4 ps, is very short. For times of the order of 1 ps or less, 
molecular motion as represented by a diffusion coefficient cannot be 
defined. Hence, in the system succinonitrile/water, the condition, 
t >> 0, required for the validity of the asymptotic formulae, is established 
rapidly* This means that the theory should be applicable as soon as the 
initial nucleation phase is complete and the surface between the droplet and 
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the continuous phase is sufficiently established for an interfacial energy to 
be defined. 


It is useful to compare the phase separation of metallic alloys 
with the system succinonitrile/water . For liquid alloys, the interfacial 
energy is about a factor of 10 higher, which increases a by a factor of 10, 
assuming room temperature. On the other hand, according to Eq. 

(3.8.5), 0 increases by a factor of 100. Nevertheless, because of Eq* 
(3.8.13), droplets of 1 ym size can be expected to be observed after 8.9 sec 
at room temperature, which is a factor of 10 sooner. On the other hand, in 
solid non-equilibrium metal alloys, the diffusion coefficient should be about 
lO" 11 cm 2 /sec, implying that 0 increases by about a factor of 10^ 2 over 
succinonitrile/water. The time to form grains of 1 ym extends to 8.9 x 10^ 
sec, implying that these non-equilibrium alloys are metastable over times 
substantially exceeding their expected useful lifetimes. 


84 


REFERENCES 


1* W. Ostwald, Z, Physik. Chem. 34, 495 (1900). 

2. F. Rosenberger, Fundamentals of Crystal Growth I , Springer-Verlag, 
Berlin, 1979, Chap. 5. 

3. I. M. Lifshitz and V. V. Slyozov, J. Phys. Chem. Solids 19 , 345 (1961). 

4. C. Wagner , Z. Electrochem. 65 , 581 (1961). 

5. J. A. Marqusee and J. Ross, J. Chem. Phys . 79 , 373 (1983). 

6. J. A. Marqusee and J. Ross, J. Chem. Phys. 80, 536 (1984). 

7. M. Smoluchowski , Z. Physik. Chem. 92 , 129 (1918). 

8. J. H. Northrop and M. L. Anson, J. Gen. Physiol. 12, 543 (1929); 

R. A. Robinson and R. H. Stokes, Electrolyte Solutions , 2nd Ed., 

Butterworths, London, 1965, Chap. 10. 

9. J. K. Baird, J. Phys. Chem. (submitted). 

10. R. F. Sekerka, C. F. Jeanfils, and R. W. Heckel, "The Moving Boundary 
Problem," in H. I. Aronson, Lectures on the Theory of Phase Transitions , 
American Institute of Mining, Metallurgical and Petroleum Engineers, New 
York, 1975, pp. 117-159. 

11. J. E. Mayer and M. G. Mayer, Statistical Mechanics , John Wiley and Sons, 
New York, 1940. 


85 


APPENDIX A - THE DELTA FUNCTION AND ITS PROPERTIES 


To 

PAUL DIRAC 

who saw that it must be true, 

LAURENT SCHWARTZ 
who proved it, 

and 

GEORGE TEMPLE 

who showed how simple it could be made 

The above is from the dedication to the book by M J. Lighthill, Fourier 
Analysis and Generalized Functions, Cambridge University Press, 1964. 

The Dirac delta function has many applications in applied mathematics, 
not the least of which is its use as a representation of a point source with 
the partial differential equations of mathematical physics. 

The delta function can be represented as a limit of a sequence of func- 
tions. Consider, the sequence, defined by (other sequences are also possible). 


6 ( x) = — T"7o exp(-n 2 x 2 ), n = 1, 2, ••• . (A. 1) 

n i / £ 

7T 

The delta function is the limit 

S(x) = lim 6 ( x) = lim —ryz exp(-n 2 x 2 ) . (A. 2) 

n I / 1 

n>«> n-H» -Jr 

In Figure A. 1 , we show how 6 n (x) looks for increasing values of n. Each 
<$ n (x) gets taller and narrower as n becomes larger. In the limit n+oo, S^(x) 
gets infinitely tall and narrow (a point source), and is called the Dirac 
delta function, <5(x). 

Although <5(x) is highly singular, it has many uses when it appears 
within integrals. For example 


+00 


/ dx6(x) 


+ 00 

f dx lim 6 (x) 
J n 

—oo n->oo 


+oo 2 2 

= lim - TJ 2 ^ ^ xe n X = 1 • 

n-v<» it -00 


(A. 3) 
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Figure A. 1 • Sequence of functions representing the Dirac delta function 

in the limit n+oo. 

Summarizing, this important property is 
+ 00 

f dx6(x) - 1* (A. 4) 

—CO 

Let f(x) be an integrable function of x* Another important property is 
+00 

j dxf(x) 6 (x - x* ) = f ( x 1 ) . (A. 5) 

—00 

We demonstrate this using Eq. (A. 1) 


88 


(A.6) 


6(x - x' ) = lim " exp[-n 2 (x - x' ) 2 ] . 

n+oo it 


To be specific, let f(x) = e x , then 


+00 


J dxe 6 ( x - x' ) = 


+0 ° 2 , .*2 

, . n f . x -n (x - x' ) 

lim J dx e e 

n+oo TJ -00 


= lim 


2,2 +oo 22 2, .. 

n -n x' t , -n x + (2n x' + 1)x 


n+oo ir 


1/2 


f dx 


2 .2 2,2n 2 x'+1^ +<» 2 r i, (2n 2 x‘ + 1) 

n -n x' n ( — ) , . -n x - '/ 2 ^ 

lln — e • e ' 2n 2 ' f 0* * ' n 2 

n>00 Jf -00 


2 .2 2,2n 2 x' + 1 % 2 1/2 

n -n x* n ) it 

lim . rr e • e ' _ 2 7 

1/2 2n n 

n->00 7 T 


e lim e 

n^oo 


(1/4n ) x' o x* 

= a — Q 


= e • e = e 


We have shown explicitly that 

+°0 , 
f dx e X fi(x - x') = e X • 

— 00 

The delta function can be generalized to three dimensions* In 
lar coordinates/ S(r - r 1 ) means 

<$(r - r’ ) = 6(x - x*)6(y - y')d(z - z')* 

The integral property expressed by Eq, (A* 5) can be generalized to 

f(r* ) = / d 3 rf (r) 6(r - r’ ) 


(A.7) 


rectangu- 


(A.8) 


(A. 9) 
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where the symbols in Eq. (A. 9) have the meanings: 


f ( r) = f(x, y, z) 
d 3 r = dxdydz . 

The integral sign in Eq. (A. 9) represents a triple integral, namely 


+00 +00 +co 

f d 3 r = f dx f dy f dz. 

— oo — oo — oo 


Combining Eqs. ( A. 8)- ( A. 1 2 ) , we can show how the right-hand side of 
is computed. 


+ 00 +oo +oo 

f d 3 rf(r) 6(r' - r' ) = f dx<5(x - x' ) J dy6(y - y' ) / dz6(z - z' 

— oo — oo — oo 


+00 +00 

= f dxS ( x -x' ) f dyS(y - y' ) f(x,y,z') 

— 00 — oo 


+oo 

= / dx6(x - x' ) f ( x,y ' ,z ' ) = f(x' ,y',z') 
— 00 


where we have used Eq. (A. 5) on each coordinate separately. 

Finally, it is worth noting that by virtue of Eqs. (A.2) and 
delta function is an even function of its argument, namely 

6 ( r - r') = 6(r' - r) . 


(A. 10) 

( A. 1 1 ) 

(A. 12) 
Eq. (A. 9) 

) f(x,y,z) 

( A. 1 3 ) 

(A.8), the 
( A. 14) 
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APPENDIX B - FICK'S FIRST AND SECOND LAWS WITH DELTA FUNCTION SINKS 


Fick's first law of diffusion is 


J(x) 


-D 


dC( x) 
dx 


9 

where C(x) is the concentration [say molecules/cm ] of the diffusing species 

2 

and J(x) [molecules/ cm sec] is the flux of that species measured across a 
plane whose normal points in the +x-direction. 

Fick's first law makes physical sense* Consider the concentration 
profiles shown in Figure B. 1 . 




(a) (b) 

Figure B. 1 . (a) Negative concentration gradient showing 

direction of molecular flux; ( b) Positive concentration 
gradient showing direction of molecular flux. 


In both cases/ J(x) points in a direction so that particle motion serves to 
smooth out the concentration gradient* 

Generalized to three dimensions/ Fick's first law reads 


J(r) = -DVC(r) 


( B* 1 ) 


where r - (x/y f z) is the position coordinate locating a point in space where 
the concentration is C(r). 
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Now, return to one dimension and consider the flux J(x) crossing a plane 
at x and another flux J(x + dx) crossing a plane a distance dx away at x + 
dx. This is depicted in Figure B.2. 



Figure B.2. Geometry for derivation of Fick*s second law. 


Inside the volume enclosed by the two planes, we allow material to be gene- 
rated by some sources (e.g., chemical reactions). For mathematical con- 
venience, we may separate these into two types: (1) Continuous sources are 

defined everywhere and are represented by continuous functions, (|>(x). The 
units of 4>(x) are, for example, molecules/cm sec. (2) Point sources are 
represented by delta functions as q5(x - x* ) where q is a constant (strength 
of the source) and x* lies between x and x + dx. If the units of x are cm, 
then the units of 6(x - x' ) must be cm so that the integral 


+00 

f dx6 ( x - x' ) = 1 
—00 


o 

is dimensionless. This implies that q has the units molecules/cm sec; hence, 
q6(x - x’) has the units molecules/cm 3 sec. 
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The distinction between a source and a sink is a matter of sign. Every- 
where <f>(x) < 0, it describes a sink. Likewise, if q < 0, q<5(x - x 1 ) is a 
point sink. 

The buildup (depletion) of material between x and x + dx can be de- 
scribed in terms of the fluxes and the sources and sinks. Let the planes at x 
and x + dx have the common cross section, A. The volume of material between 
the planes is Adx, whereas the amount of material is C(x)Adx. If we allow all 
functions to depend also upon the time t, then the net rate of buildup of 
material within this volume is (9C(x, t) /at) Adx and 


9C(x, t) 
9t 


Adx = J(x, 


t) A - J(x + dx. 


t) A + 4>(x)Adx + q<$ ( x - x' )Adx 


( B. 2 ) 


where we have allowed for the presence of both continuous and point sources. 
We expand J(x + dx, t) in a Taylor series about x, as 


J(x + dx, t) = J(x, t) + — — dx 

dX 


(B.3) 


keeping only terms up to the first order in dx. We now substitute Eq. (B.3) 
into Eq. (B.2) to obtain 


9C(x, t) , „ ,, 9 J(x, t) 

— Adx = J(x, t) Adx - J(x, t) Adx - Adx 

at ax 


+ <|)(x)Adx + q6(x - x 1 ) Adx • 


(B . 4) 


After cancellation of terms of opposite sign, common factors, and algebraic 
r e ar r angemen t 


3C( x, t) + 3J( x, t) 
3t 3x 


<J>(x) 


+ q<5(x - x' ) 


(B.5) 


Eq, (B.5) can be generalized to three dimensions: 


9 3 


0 2 . 


(B.6) 


+ v 
3t 


J(r, t) = <J> (r) + q 6 (r - r' ) 


where we have used the definition of 6 (r - r' ) given by Eq. (A. 8 ). We substi- 
tute Eq. (B.2) and obtain 


- D7 2 C(r, 


t) = 4 , (r) + qS(r - r' ) 


(B.7) 


which expresses Fick's second law. 

In the steady state/ the concentration field/ C(r, t) is no longer a 
function of time and 8C/3t = 0/ so that Eq. (B.7) becomes 

DV 2 C(r) = - 4 >(r) - q 6 ( r - r’ ) . (B. 8 ) 
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APPENDIX C - THE GREEN'S FUNCTION METHOD 


C. 1 General Principles 

The Green's function method is a technique for integrating inhomogeneous 
partial differential equations. Essentially r it is a generalization of the 
method of variation of parameters used to solve inhomogeneous ordinary dif- 
ferential equations. 

In steady state diffusion, we are faced with the solution of a partial 
differential equation of the form 

DV 2 C(r) = f ( r) (C.1) 

where f(r) is a source term which, by virtue of its presence, makes Eq. (C.1) 
inhomogeneous. There must also be some conditions specifying the value 
of C(r) on the boundaries. 

C. 2 Green's Problem 

In order to solve Eq. (C.1) by the Green's function method, we suppose 
that we may solve exactly the equation 

DV 2 G(r|r* ) - 6(r - r* ) (C.2) 

where G(r|r' ) is called the Green's function, whose values we are free to 
choose on the boundaries. Because of the singular properties of the delta 
function, Eq. (C.2) is homogeneous everywhere except at the point r = r' • Eq. 
(C.2) is known as Green's problem. 

We proceed to solve Eq. (C.2) for G(r|r'). To accomplish this, we make 
use of the divergence theorem, which states 

/ d 3 r V . F(r) = / dS • F(r) (C.3) 

where F( r) is a vector-valued function of the position r and where in rectang- 
ular coordinates, r - (x,y,z) . The triple integral covers the entire volume 
enclosed by the boundary S. 
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Note that 


v^ 2 G(s|r') = V ? • (7 G(5|r , )l. 

Integrate Eq. (C.6) over all space starting from the origin at O'. 


(C.7) 


/ <S 3 56 (£) = 1 = f d 3 £ D • fv ? G(5|r' )) 

= D / dS ? • V ? G(ejr') (C * 8) 

where we have used Eq. (C.3) with F(r) = G(g|r* )• We let our volume of 

interest be a sphere centered on 0 1 • The radius of this sphere is g = le! 

2 

which has surface area 4 tt£ • The normal to the surface of this sphere points 
in the ^-direction. We intend for 6(g) , to represent a point source (iso- 
tropic) centered on 0 1 ; hence, G(g|r* ) must also be isotropic (i.e., depending 
only upon g = Ul and independent of angles). The surface integral in Eq. 
(C.8) is then equal to 

D / dS • V G( c|r' ) = 4ir£ 2 D- 8 - - ( f [— ■■ <C.9) 


where 3G/3g is the gradient of G in the g-direction (the radial direction for 
spherical polar coordinates centered on O')* Putting Eq. (C.9) into Eq. 
(C.8), we obtain 


8G(£lr' ) 
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1 

4ttD£ 2 


Eq. (C.10) can be integrated to read 


(C. 10) 


+ c (c -”> 

where C is a constant. We take as boundary condition, G(oo|r* ) = 0, from which 
we conclude 


G(r|r" ) 


1 

4ttD| r - r' | 


( C. 12) 


where we have used Eq. (C.4) and g = Ul = I r " r * | • The gradient of 
G( r | r' ) is 
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V G(r r’ ) = 


_L [i i- + j L + i i_i 

4ttD ' 3x J 3y 3z 1 


f(x-x' ) 2 + (y-y') 2 + (z-z') 2 ] 


1 r (x-x') i + (y-y 1 ) j + (z-z') k ■ 

4ttD 3/2 

((x-x' ) + (y-y') + (z-z') ] 


(C.13) 


which, in the limit that x, y, z all go to infinity, is zero. Hence, the 
Green's function expressed by Eq. (C.12) satisfies the boundary conditions 


lim G(r|r' ) = lim V G(r|r' ) = 0. 


(C.14) 


Note further that because of the absolute value sign in Eq. (C.12), G(r|r' ) is 
symmetric upon interchange of r and r' , namely 


G(r|r' ) = G(r' |r) . 


( C. 15) 


This is called the reciprocity relation for the Green's function and is due to 

2 

the fact that the operator V is self-adjoint. 


C. 3 Green's Theorem 

We propose to solve Eq. (C.1) by expressing C(r) in terms of some inte- 
grals involving G(r|r'), f(r), and the values of C(r) specified on the bound- 
aries (boundary conditions). To begin, we reexpress Eqs. (C.1) and (C.2) as 
partial differential equations in terms of the variable, r* , 


DV' C(r' ) = f ( r' ) 


(C. 16) 


DV' 2 G(r' | r) = 6(r' - r) 


( C. 17) 


where V' is the Laplacian operating on the variable r' • Now multiply Eq. 
(C.16) through by G(r'|r) and Eq. (C.17) through by C(r' ) and subtract. The 


result is 


G(r'|r) DV ,2 C( r' ) - C(r' )DV ,2 G(r|r' ) = G(r'|r)f(r') - C(r ')6(r'- r) . (C.18) 
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The left-hand side of Eq. (C.18) may be simplified by using the vector identi- 
ty 


V 1 • (A ?'B) 


V'A.V'B + A V'B 


(Cm 19 ) 


where A(r* ) and B(r* ) are functions of r* • Eq. (C.19) can be rearranged to 
read 


A V ,2 B = -V'A.V'B + V' »(AV'B) . (C.20) 

Using Eq. (C.20) with A = G and B = C, the first term on the right-hand side 
of Eq. (C.18) may be written 

G(r' |r)DV 2 C(r' ) = DV' • (G(r' |r) V'C(r' ) ) - DV'G(r' | r) . V'C( r' ) (C.21) 

and using Eq. (C.20) with A = C and B = G, the second term on the left-hand 
side may be written 

C(r’ )DV' 2 G(r’ |r) = DV' »(C(r' )V'G(r' |r) ) - DV'C(r' ) •V'G(r’ |r) . (C.22) 

Substituting Eqs. (C.21) and (C.22) into Eq. (C.18), we obtain 

DV , *[G(r , |r)V , C(r') - C( r' ) V'G(r' |r) ] = G(r' |r)f (r' ) -C(r')6(r' - r) . (C.23) 

We next integrate both sides of Eq. (C.23) over the volume of all space, 
3 

i.e., / dr', and apply the divergence theorem to he left-hand side. The 
result is 


D f dS* • [G(r' |r)V'C(r' ) - C(r' )7'G(r' |r) ] 


= / d 3 r'G(r’ jr)f(r' ) - / d 3 r’C(r' ) 6(r' - r) 


( C. 24) 
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On the basis of Eq. (A. 9), the second term on the right-hand side of Eq. (C.8) 
( note that the roles of r and r* have been interchanged) is just equal to 
C(r). Hence, Eq. (C.24) may be rewritten 


i 


C(r) = D / dS'. [C( r' ) V'G(r' | r) - G(r' |r) V'C(r’ )] 


+ / d 3 r'G(r’ |r)f(r') 


If we let the surface. S', recede infinitely far away, i.e., we 

may employ Eq. (C.14) to evaluate the surface integral in Eq. (C.25). We 
conclude that so computed this integral is zero. Hence, we obtain finally 
frcm Eq. (C.25) the result 


C(r) = f d 3 r , G(r|r* )f (r 1 ) 


where we have used the reciprocity of G(r|r f ) expressed by Eq. (C.15). Pro- 


vided that the integral on the right-hand side of Eq. (C.26) may be computed. 


this serves to solve Eq. (C. 1) 


APPENDIX D - FOURIER TRANSFORMS 


D. 1 General Principles 

Let <(>(x) be a function of the variable x where -<» < x < + <», The 

A 

Fourier transform/ <J>(k) of <|>(x) is 


<(>(k) 


/ dx d (|>(x) 


(D.1) 


where k is the transform variable and i = (-1) 1 / 2 . Starting with Eq. (D.1)/ 

* 

Fourier showed that <f>(k) and $(x) are also related by 


[x) = L I dk e ikx <f>(k) . 


(D.2) 


Equation (D.1) may be generalized to three dimensions. Consider 


<(>(k.j , k 2 , k 3 ) = / dx 


+» -ik^x +oo “ik 2 y +oo -ik 3 z 

e f dy e f dz e 4>( x,y,z ) (D.3) 


or 


-ik«r 


4 > ( lc) = / d r d 4 >(r) 


(D.4) 


where the triplet (k^, k 2 , k 3 ) is treated as a vector, k = (k 1 , k 2 , k 3 ) and r 
= (x, y, z) . Equation (D.2) generalized to three dimensions becomes 


+ GO 


1 ■ - ik x 1 +oo ik y +oo ik z A 

4»(x,y,z) = (— ) / dk 1 e (—1 / dk 2 e (— 1 f dk 3 e <(>(k) (D.5) 
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or 


4> (r) = (•— ) / d 3 k d lk * r J(k) 


(D.6) 


D. 2 Fourier Integral Representation of the Delta Function 
In Eq. (D.4), let 4> { r ) = <5(r - r' ) , then 


7,,, c , 3 -ik»r -ik«r' 

<f> ( k) = J d r e 6(r - r' ) = e 


(D.7) 


We substitute Eq. (D.7) into Eq. (D.6) to obtain 


SOr-r*)^] / d 3 kd- ik * (r - r '> 


(D.8) 


Equation (D.8) is the very valuable Fourier integral representation of the 
delta function* It is much more useful than the representation expressed by 
Eq. (A. 2)* 

An integral representation of the delta function in k-space, namely 
6(k - k' ) , is obtained if in Eq. (D.6) we let <}> ( k) = 6(k - k* ). We find 

♦ (r) = / d 3 k e lk * r 6(k - k') = (-t-1 e lk ' * r . (D.9) 


Now substitute Eq. (D.9) into Eq. (D.4) 


<5(k - k’ ) 


- (5) / • 


-i(k - k' ) .r 


(D . 10) 


Equation (D.10) is an integral representation of 6(k - k* ) • 

The delta function in Fourier integral space has properties analagous to 
those it exhibits in real space, namely 
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(D.11) 


/ d 3 k’ 6(k' - k) = 1 
/ d 3 k* 6(k' - k)f(k' ) = f(k). (D. 12) 

D. 3 Convolution Theorem of the Fourier Transform 

The function f(x) is said to be the convolution of the functions g( x) 
and h(x) if 

+ 00 

f ( x) = / dx' g(x-x* )h(x’ ) = g*h( x) (D.13) 


where g*h(x) is the short hand notation for the convolution of g(x) and h( x) 
and is defined by the integral to the left. 

The Fourier transform of f(x) is 


A +°0 . +00 

f(k)= J dx e 1 X f dx' g( x - x' ) h( x' ) 


In the integrand, let £ = x - x' , d£ = dx, so that 


f (k) 


/ d5 / dx' e lk( ^ + X ^g(5)h(x' ) 


- / de e‘ ik5 g(C) 


+ OO 

/ dx' 


— CO 


-ikx' 

e 


h( x' ) 


f(k) = g(k) h(k) 


where we have used 


(D. 14) 


(D. 15) 


g(k) 


f d? e" lkC g(C) 


(D.16) 
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h(x’) 


(D.17) 


h(k) 


+ 00 

= f dx' 


-ilex' 

e 


— oo 


The above results may be generalized to three dimensions , namely 
f ( r) = f d 3 r' g(r - r' )h(r' ) = g*h(r) 


( D* 1 8 ) 


and 

f(k) = g(k) h(k) ( D . 1 9 ) 

where the transforms are defined in the sense of Eq. (D.4). We conclude that 
the transform of a convolution is equal to the product of the individual 
transforms. 

2 

D.4 Fourier Transform of V C(r) 

Let C(r) be a scalar-valued function of the coordinates r = (x, y, z) . 
The Laplacian of C(r) expressed in Cartesian coordinates is 


V 2 C(r) 


2 2 3 

8 C a_c 3_C 

2 2 2 
3x 3y 3z 


( D. 20 ) 


2 2 

For simplicity/ we compute explicitly only the transform of 3 C/3x . We 
integrate by parts. 


2 

r 3 -ik«r 3 C 
fdre — 

dx 


+oo 

f dx 


— oo 


"ik^x +«> ~ ik 2 y +00 " ik 3 Z 9 2 c 

e f dy e f dz e — j 

-oo -00 3x 


+oo -ik y +oo -ik z +oo -ik.-x 2 

= f dy e / dz e f dx e — - (D.21) 

— OO —00 —00 3 X 
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+<» -ik y +oo -ik z -ik x +«> -ik„x +oo 


f dy e f dz e ||e 

-oo “oo 


«* ’ S ♦*»,[' 'Cl 


l 


+oo -ik„x 


- kj / dx e 1 C} 


The boundary conditions are often 


lira C(x, y, z) = lira = 0 (D.22) 

X->±oo X+±co 


in which case Eq. (D.21) simplifies to 


e _3 -ik*r 3 2 C 
/ d r e — o 

ax 



+00 -ik X +00 

/ dx e f dy 



f dz 


— 00 


-ik 3 y 

e C(x,y,z) 


(D.23) 


or 


, ,3 -ik.r a 2 C 
/ d r e — 

3x 


2 * 

k 1 C(k) 


( D. 24 ) 


where 


C(k) = / d 3 r e ik>r c(r) . 


(D.25) 


Generalizing to three dimensions, the transform of Eq. (D.20) is 
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(D.26) 


f j 3 -ik«r 2 . . , 2 * . 

J d r e V C(r) = -k C(k) 


where 


k 2 = k 2 + k 2 + k 2 = |k | 2 . (D.27) 

D. 5 Fourier Transform of G(r) = -(4iTDr) ^ 

In Appendix C.2, we showed that the solution to Green’s problem for the 
steady state diffusion operator was given by Eq. (C.12). If in Eq. (C.12) 
r' = 0, we obtain 


G(r) 


1 

4irDr 


(D.28) 


It is useful to have the Fourier transform of Eq. (D.28). To compute this, we 
substitute Eq. (D.28) into Eq. (D.4). The result is 


G(k) = -(1/4irD) / d 3 r e" lk * r (1/r) . 


(D.29) 


We evaluate the integral using spherical polar coordinates. Let k and the 
polar axis (z-axis) point in the same direction (see Figure D.1). We obtain 
for the inner product, k«r, 

k*r = kr cos0, (D.30) 

and for the volume element 
3 2 

d r = dr r d0 sin0 d<J> • (D.31) 


106 


Figure D. 1 


Geometry for computing the integral in Eq. (D.29) 


Transformed to spherical polar coordinates, Eq. (D.29) becomes 


TT 


G(k) = -(1/4 ttD) f dr r (1/r) f d0 sin0 e 


-ikr cos0 


2tt 

f <4 . 

O 


(D. 32) 


If we let x = cos0 , and dx = -sin0 d0 


G(k) = — ( 1/4nD) / dr r 2 (1/r) 
o 



dx e 


-ikrx 


2 TT 

f d<j) 


o 


= -(1/4ttD) f dr r 2 ( 1/r) (2/kr) (sin kr)(2ir) 
o 


00 

(1/kD) f dr sin kr . 


(D. 33) 



The remaining integral in Eq. (D.33) may be evaluated by writing sin kr in 
terms of exponentials and inserting an integrating factor as 


00 00 ikr -ikr 

/ dr sin kr = / dr (- ^ ) 


= (-t-) {lim f dr e lkr ‘ er - lim / dr 

" £>0 O £->“0 O 


-ikr-sri 
e } 


ikr-er « -ikr-er « 

= film r — tt 1 - lim r e ,.. - - , 1 } 

^2i 1 ^ 1 ikr-e 1 „ -(ikr+e) ' 

e+o o e-»0 o 


t 11 " (ifel * I - 'A . 

£->■0 £>0 


We substitute Eq. (D.34) into Eq. (D.33) and obtain our desired result 


G(k) = - 1/k D . 


(D. 35) 
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APPENDIX E - SPECIFICATION OF EXPONENTS 


E.1 Proof of y-j = x 

We first consider the limit of W(Z, t) for large t. To do this, we 
substitute Eq. (3.3.6) into Eq. (3.3.22) and write 


W(Z,x) 


= t" X (1/Z) 


-Yi 


[0 1 T 


t" X (1/z)] [1 + eZ] . 


(E.1. 1) 


The functional form of W(Z, x) for large x depends upon the magnitude of as 
compared to x. Letting x-k» in Eq* (E*1*1), we find the results 


W(Z /T ) 


-( 1/z) 2 (1+eZ) t~ 2x 

(1/z)(o 1 - ( 1/z) ) ( 1+eZ) x 
-(x+y ) 

(1/z)a^x 


if y 1 > x 

if y 1 = x 


if Y 1 < x 


(E.1.2) 


Equation (E.1.2) in turn may be written as a single equation 
w(z, T ) = t - t1 V(Z) 


(E. 1.3) 


where both n and V(Z) depend upon the 
Substitute Eq* (E.1*3) into Eq. 


relative magnitudes of y^ and x. 
(3.3*19) and obtain 




} . t -(S^) 31F o (Z)V|Z)1 


az 


(E.1.4) 


Equating exponents of x in Eq. (E.1.4) , we find 

T) = 1 “X • (E. 1.5) 

We next consider the time derivative of the super saturation a(x) given 
by Eq. (3.2.9), which we reproduce below 
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“ 3 

o(t) = CT(0)[1 - / da F(a, x )a ] . 
o 


(E.1.6) 


The derivative of Eq. (E.1.6) with respect to x is 


3cr(x) 

3t 


/ ^ % r - 3F ( a / T ) _ 3 

- a(0) J da — — a 


(E.1.7) 


We may write Eq. (3.2.18) as 


3F( a ,x) = 
3x 


3 [F(a,x)W(a,x)1 
3a 


(E. 1.8) 


and substitute into Eq. (E.1.7) to obtain 


3a(T) 

3t 


ff(0) / 

o 


da a 


3 3 [F( a, t)W( a, T ) ] 

3a 


(E.1.9) 


Using F(», T ) = 0, Eq. (E.1.9) may be integrated by parts. The result is 


3q( t) 
3t 


-3a(0) 


2 

/ da a F( a,x)W( a,x) • 
o 


(E.1.10) 


Using Eq. (3.3.2) to convert Eq. (E.1.10) to an integral over the canonical 
length Z, we obtain 


"fx” = ~ 3<7 ^ j3X / dz z ^F(Z,x)W(Z,x) • 


(E. 1.11) 


In the limit T+°°, we may replace W(Z,x) by Eq. (E.1.3) and F(Z,x) by Eq. 
(3.3.8). In this limit Eq. (E.1.11) becomes 
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— • — = -3a(0) T _(x+n) / dZ z 2 F (Z)V(Z) . 

9t : o 


(E. 1.12) 


On the other hand, differentiating Eg. (3.3.6), we obtain 


? 


3a(t) 

3x 


-(y^D 

- y i a i T 


(E.1.13) 


t 


and substituting Eq. (E.1.13) in the left-hand side of Eq. (E.1.12) while 
using Eq. (E.1.5), we have 


3a(0) / dZ ZF (Z)V(Z) = y.cr.T 
J o II 


-yi 


(E.1.14) 


~ y i 

Since y 1 > 0, x *0 as I* 00 * The left-hand side of Eq. (E.1.14), however, 
independent of t* If the left-hand side is to be zero, the integrand must 
have compensating positive and negative portions. In as much as F q (Z) is a 
distribution function of grain sizes, it is positive definite, which implies 
that V(Z) must possess both positive and negative swings. Referring to Eq. 
(E.1.2), we find that this is true only for 


is 


v(z) = ( 1/z) ( cr 1 - (1/z))<1 + e z) 


(E.1.15) 


which imples 

y^=x. (E. 1 • 16) 

E.2 Proof of x = 1/3 

Since y-j = x, we find by comparing Eqs. (E.1.2), (E.1.3), and (E.1.16) 

that 
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n = 2x 


(E.2.1) 


Substituting Eq. (E.2.1) into Eq. (E.1.5), we conclude 


x = 1/3 . 


(E. 2.2) 


t 


I 


/ 
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APPENDIX F - EXPANSION OF 1/w(Z) IN PARTIAL FRACTIONS 
Rewriting Eq. (3.5.8), we may use partial fractions. 


1 


Al , A 2 , B 1 


3w(Z) (Z-Z ) 2 (Z+3/Z 2 ) (Z-Z ) 2 (Z-Z ) (Z+3/Z 2 ) 

o o o o o 


The coefficients A^ , A 2 , and are obtained by evaluating 


2 2 

(Z-Z ) Z 

», - 1- — 2 °— — ) z . 


(Z-Z ) 2 (Z+3/Z 2 ) Z Z o (Z 3 +3) 
o o o 


2 2 3 3 

. (Z-Z ) Z Z J (Z J +6) 

A = r 2 l = ° 2 

2 dZ ,rr rr X 2 /„ 2 \ J Z=Z . „3 J _- , 2 

(Z-Z ) (Z+3/Z ) o (Z +3) 

O O O 


(Z+3/Z 2 ) Z 2 

l 1 = I 2 1 _ 2 

(Z-Z ) 2 (Z+3/Z 2 ) Z - 3/Z 0 
o o 


3 2 

(Z +2T 

o 


respectively. Substitution of Eqs. (F.2)-(F.4) yields 


3Z 


1 


W(Z) (Z 3 +3) (Z-Z.) 2 

o o 


3z 3 (z 3 46) 
o o 

3 2 

(Z +3) 

O 


27 


(Z-Z ) 
o 


3 2 

(Z>3) 

o 


(Z+3/Z ) 
o 


(F.1) 


(F.2) 


(F.3) 


(F.4) 


• (F.5) 
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16. Abstract 

When a two-component system is cooled below the minimum temperature for its stability, it separates into two or more immiscible 
phases. The initial nudeation produces grains (if solid) or droplets (if liquid) of one of the phases dispersed in the other. The dynamics by 
which these nuclei proceed toward equilibrium is called Ostwald ripening. We shall assume that the nuclei are spheres, thereby treating 
the case of droplets rigorously and the case of grains approximately. 

The dynamics of growth of the droplets depends upon the following factors: (1) The solubility of the droplet depends upon its radius and 
the interfacial energy between it and the surrounding (continuous) phase. The equation governing this phenomenon is associated with the 
names of Gibbs and Kelvin, and we derive it for the case of two components. There is a critical radius determined by the supersaturation 
in the continuous phase. Droplets with radii smaller than critical dissolve, while droplets with radii larger grow. (2) The droplets concen- 
trate one component and reject the other. The rate at which this occurs is assumed to be determined by the interdiffusion of the two 
components in the continuous phase. Associated with this diffusion is advection, which is accounted for in our theory to lowest order in 
the center of volume velocity of the two components. Moreover, when the droplets occupy a finite fraction of the total volume of the 
system, there occur effects of inter-droplet competition on the concentration field in the continous phase. This increases the overall rate 
of growth of the droplet phase. Using diagrammatic techniques and rigorous statistical mechanics, we show that this effect depends upon 
the square root of the volume fraction associated with the grains, and that it accelerates both the rate of growth and the rate of dissolu- 
tion of droplets. (3) The Ostwald ripening is constrained by conservation of mass; e.g., the amount of material in the droplet phase plus the 
remaining supersaturation in the continuous phase must equal the supersaturation available at the start. (4) There is a distribution of 
droplet sizes associated with a mean droplet radius, which grows continuously with time. This distribution function satisifes a continuity 
equation, which we solve asymptotically by a similarity transformation method. 

The final results of our theoretical analysis consist of the following quantities determined as functions of time: (1) the supersaturation 
remaining in the continous phase, (2) the total number of droplets in the continuous phase, (3) the distribution of droplet sizes, (4) the 
critical droplet radius, (5) the mean droplet radius, and (6) the maximum droplet radius. Each of these is also determined as a function of 
the volume fraction occupied by the droplets. We apply these results to the system, succinonitrile/water, which separates into two immis- 
cible liquid phases below its consolute temperature at 58 °C. This system is particularly attractive because it is transparent to visible 
light, so the rate of formation of the droplet phase can be followed photographically. In general, the droplet phase has a different density 
than the continuous phase. Separation of the two by sedimentation competes with the process of Ostwald ripening. In the acceleration- 
free environment of an Earth orbiting laboratory, the Ostwald ripening in the system succinonitrile/water should be observable without 
interference by sedimentation or any other effects associated with gravitation. We provide a complete set of specifications for the design 
of such an experiment- 
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